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ABSTRACT 



We present a detailed analysis of the gas conditions in the H2 luminous radio galaxy 3C 326 N at z~0. 1 , which has a low star-formation 
rate (SFR~0.07 M ;;j yr ) in spite of a gas surface density similar to those in starburst galaxies. Its star-formation efficiency is likely 
a factor ~ 10-50 lower than those of ordinary star-forming galaxies. Combining new IRAM CO emission-line interferometry with 
existing Spitzer mid-infrared spectroscopy, we find that the luminosity ratio of CO and pure rotational H2 line emission is factors 
10-100 lower than what is usually found. This may suggest that most of the molecular gas is warm. The Na D absorption-line profile 
of 3C 326 N in the optical suggests an outflow with a terminal velocity of ~ — 1 800 km s~ 1 and a mass outflow rate of 30-40 Mq yr~ , 
which cannot be explained by star formation. The mechanical power implied by the wind, of order 10 43 erg s _1 , is comparable to the 
bolometric luminosity of the emission lines of ionized and molecular gas. To explain these observations, we propose a scenario where 
a small fraction of the mechanical energy of the radio jet is deposited in the interstellar medium of 3C 326 N, which powers the outflow, 
and the line emission through a mass, momentum and energy exchange between the different gas phases of the ISM. Dissipation times 
are of order 10 7-8 yrs, similar or greater than the typical jet lifetime. Small ratios of CO and PAH surface brightnesses in another 
7 H2 luminous radio galaxies suggest that a similar form of AGN feedback could be lowering star-formation efficiencies in these 
galaxies in a similar way. The local demographics of radio-loud AGN suggests that secular gas cooling in massive early-type galaxies 
of > 10 Mq could generally be regulated through a fundamentally similar form of 'maintenance-phase' AGN feedback. 

Key words. Galaxies - ... ... 



1. Introduction 

Molecular gas plays a critical role for our growing understand- 
ing of galaxy evolution. It often dominates the mass budget of 
the interstellar medium in galaxies, and is most closely related to 
the intensity at which galaxies form stars (e.g., |Kennicutt|1998] l. 
Being strongly dissipative, it is also particularly susceptible to 
the astrophysical processes that drive galaxy evolution - inter- 
actions, or feedback from starbursts and AGN - and therefore 
plays a key role for our understanding of how these processes 
regulate star formation and galaxy assembly. 

It has only recently been recognized that powerful AGN 
may play a significant role in regulating galaxy growth over 
cosmological timescales by suppressing gas accretion and star 
formation (e.g., |Silk & Rees||1998||Friaca & Terlevich|[ l998 
Scannapieco & Oh|2004[|Springel etal.|2005| [Bower et al.|2006 
Croton et al.||2006l |Ciotti & Ostriker||2007| ^Merloni & Heinz| 
2008) >. Such AGN 'feedback' would help resolve some of the 
remaining discrepancies between hierarchical models of galaxy 
evolution - implying a rather gradual assembly of massive galax- 
ies - and observations, which suggest that massive galaxies 



formed most of their stars at high redshift, whereas star forma- 
tion at later epochs was strongly suppressed. Observationally, a 
picture is emerging where radio jets may play a large role in 
transforming the energy ejected by the AGN into kinetic and 
thermal energy of the interstellar medium of the host galaxy. 
Observations of radio-loud AGN (e.g., |Heckman et al.|1991b|a 



Morganti et al.||2005] |Emonts et al.||2005| |Best et aO 2005 



2006; Nesvad baeT al.i2006| |2007l |McNamara & Nulsen|2007 
Nesvadba Tt al ||2008||Holt et al.||2008||Baldi & Capettifl2008 
|Fu & Stockton|2009||Humphrey et al.|2009)> and a large number 



of hydrodynamical simulations (e.g., Krause 2005; Saxton et al 
12005] IHeinz et al.||2006l ISutherland & Bicknell||2007l Merloni 



|& Heinz||2007[ |Antonuccio-Delog u & Silk 2008 ) suggest that 
radio-loud AGN inject a few percent of their mechanical energy 
into the ambient gas, parts of which produce significant outflows 
of warm gas (|Morganti et al ||2003| [2005] |Emonts et al.|[200T 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer 



Nesvadba et a l.||2006[ |20071 |Holt et al.||2008l |Nesvadba et al.| 
2008 ; Fu & Stockton 2009 ). However, most previous studies fo- 
cused on the warm and hot gas at temperatures > 10 4 K, and did 
not address the impact on the molecular phase, which is a seri- 
ous limitation if we want to understand how the radio-loud AGN 
may regulate star formation in the host galaxy. 
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Observations with the Spitzer IRS spectrograph recently re- 
vealed a significant number of "H2-luminous" galaxies, where 
the molecular gas does not appear associated with star formation 
([Appleton et d.|2006||Egami et al.|2006l|Ogle et al.|2007|[2008l 
2010, de Messieres et al. 2009; Sivanandam et al. 2009, see also 
Haas et al.||2005a|n The mid-infrared spectra of H2-luminous 



galaxies are dominated by bright, pure rotational emission lines 
of warm molecular hydrogen, (J£{Hq) — 10 40 - 10 43 ergs~'), 
while classical star-formation indicators like a bright infrared 
continuum, mid-infrared lines of [Nell] and [Nelll], and PAH 



bands are weak or absent. Interestingly, Ogle et al. (2010 1 find 



that 30% of their radio-loud AGN taken from the 3CR are H 2 lu- 
minous, suggesting this may be a common phenomenon which 
could be related to interactions with the radio source. 

To test this hypothesis and to evaluate possible consequences 
for the gas properties and star formation in the host galaxy, we 
have started CO emission-line observations of H2 luminous ra- 
dio galaxies with the IRAM Plateau de Bure Interferometer. Our 
goal is to constrain the physical properties and masses of the 
multiphase warm and cold gas in these galaxies, and to mea- 
sure the gas kinematics. Here we present a detailed analysis of 
the multiphase gas content, energetics, and dissipation times of 
the H 2 -luminous radio galaxy 3C 326 N at z=0.09 (Ogle et al 
2007| [20081 12010| , which has particularly high H 2 /PAH ratios 
This analysis is based on our new CO(1-0) observations, as well 
as existing mid-infrared Spitzer and SDSS optical spectroscopy. 
Specifically we address three questions: What powers the H 2 
emission in this galaxy? What is the physical state of the molec- 
ular gas, and perhaps most importantly, why is 3C 326 N not 
forming stars? 

We find that the interstellar medium of 3C 326 N has very 
unusual physical properties, where the warm molecular gas may 
dominate the overall molecular gas budget (<Q and where the 
emission-line diagnostics suggest that the molecular as well as 
the ionized gas may be mainly excited by shocks (f|3]l giving rise 
to luminous line emission at UV to mid-infrared wavelengths. 
We also identify a significant outflow of neutral gas from Na D 
absorption profiles, which cannot be explained by star formation 
(Q. We propose a physical framework in which these obser- 
vations can be understood as a natural consequence of the en- 
ergy and momentum coupling between the gas phases, which is 
driven by the mechanical energy injection of the radio jet (jj6). 
This scenario is an extension of the classical 'cocoon' model 
(e.g., |Scheuer|1974||Begelman & Cioffi|1989] ) explicitly taking 
into account the multiphase character of the gas, with an em- 
phasis on the molecular gas. We use our observational results 
to quantify some of the parameters of this scenario, including, 
perhaps most importantly, the dissipation time of the turbulent 
kinetic energy of the gas, and find that it is self-consistent and in 
agreement with the general characteristics of radio-loud AGN. 

3C 326 N shows evidence for a low star-formation efficiency 
leading to a significant offset from the Schmidt-Kennicutt re- 
lationship of ordinary star-forming galaxies by factors 10—50, 
which is similar to other H2 luminous radio galaxies with CO 
observations in the literature, as would be expected if 3C 326 N 
was a particularly clear-cut example of a common, underlying 
physical mechanism that is throttling star formation (|7}. Long 
dissipation times suggest that the gas may remain turbulent over 
timescales of 10 7-8 years, of order of the lifetime of the radio 
source, or perhaps even longer ({ 6.3 1, while the energy supplied 



by the radio source may be sufficient to keep much of the gas 
warm for a Hubble time (S|8]) as required during the maintenance 
phase of AGN feedback, assuming typical duty cycles of order 
10 8 yrs. 

Throughout the paper we adopt a Ho =70 km s -1 , £Im = 0.3, 
£Ia =0.7 cosmology. In this cosmology, the luminosity distance 
to 3C 326 N is D| 2av = 411 Mpc. One arcsecond corresponds to 
a projected distance of 1.6 kpc. 



2. Observations 

2. 1. The H 2 luminous radio galaxy 3C 326 N 

We present an analysis of the powerful H2 luminous radio galaxy 
3C 326 N at z~0.1 ( |Ogle et al.|2007||20T0] >. Pure rotational, mid- 
infrared H2 lines in 3C 326N have an extraordinary luminosity 
and equivalent width ( |Ogle et al.|2 007). The total emission-line 
luminosity of 3C 326N is L(H2) = 8.0 ±0.4 x lO^ergs" 1 (in- 
tegrated over the lines S(l) to S(7)), corresponding to 17 ±2% 
of the infrared luminosity integrated from 8-70 fim. This is 
the most extreme ratio found with Spitzer so far. The H2 line 
emission is not spatially or spectrally resolved, implying a size 
<4" (~6 kpc at the distance of the source), and a line width 
FWHM^SOOkms- 1 . 

3C 326 N is remarkable in that it does not show evidence for 
strong star for mation, in spite of a sign ificant molecular gas mass 
of - 10 9 M Q ( jOgleetal.|2007l|2010| ). This mass corresponds to 
a mean surface density of warm H 2 within the slit of the short- 
wavelength spectrograph of Spitzer of 3OM pc~ 2 , a few times 
larger than the total molecular gas surface density of the molec- 



ular ring of the Milky Way (Bronfma n et al.|fl988| >. Generally, 
galaxies with similarly large molecular gas masses show vigor- 
ous starburst or AGN activity. However, in 3C 326 N the lumi- 
nosity of the PAH bands and 24 ftm dust continuum suggest a 
star-formation rate of only ~ 0.07 M Q yr _1 , about 2% of that 
in the Milky Way. Given this large amount of molecular gas and 
the high mass surface density, the bolometric AGN luminosity 
of 3C 326 N is also remarkably low. 

3C 326 N has a Mpc-sized FRII radio source, which is rela- 
tively weak for this class, with a radio power of 2.5 x 10 26 W 
Hz at 327 MHz. There is some confusion in the literature 
regarding which galaxy, 3C 326 N or the nearby 3C 326 S is 
associated with the radio lobes. Both candidates have detected 
radio cores and the extended radio lobes make it difficult to 
uniquely associate one or the other component with the radio 
source. Rawlings et al. ( 1990| l argue that 3C 326 N is the more 
plausible candidate, having the brighter stellar continuum (and 
greater stellar mass, |Ogle et al.||2W7] ). Only 3C 326 N has lu- 
minous [OIII]A5007 line emission consistent with the overall 
relationship between [OIII]A5007 luminosity and radio power 
([Rawlings et al .|1990 1. The c ore of 3C 326 N appear s unresolved 
at 8.5 GHz with a 2 beam (Rawlings et al.(l990 i. We discuss 
the age and kinetic power of the radio source of 3C 326 { 4.2 



Ogle et al. |2007 2010 1 propose to introduce a new empirical 
classification based on this H2 excess, and refer to such targets as 
"Molecular Hydrogen Emission Galaxies", MOHEGs 



2.2. CO millimeter interferometry 

CO(1-0) emission-line observations of 3C 326 were carried out 
with the IRAM Plateau de Bure Interferometer (PdBI) in two 
runs in January and July/August 2008 with different configura- 
tions. In January, we used the narrow-band and dual-polarisation 
mode, corresponding to a band width of 950 MHz, or 2740 km 
s , with a channel spacing of 2.5 MHz. The 6 antennae of the 
PdBI were in the BC configuration at a central frequency of 
105.85 GHz, corresponding to the observed wavelength of the 
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CO(l-O) line at a redshift of z=0.090. The final on-source in- 
tegration time for this run was 6.7 hrs (discarding scans with 
atmospheric phase instabilities). The FWHP of the synthesized 
beam in the restored map is 2.5" x 2.1" with a position angle of 
23°. Only data from this run were used for the CO emission-line 
measurements. 

The presence of a 3 mm non-thermal continuum from the ra- 
dio source and relatively narrow receiver bandwidth made it dif- 
ficult to rule out a possible contribution of a very broad CO(l-O) 
line (FWHM>1000 km s" 1 ) from the January 2008 data alone. 
We therefore re-observed 3C 326 N in July and August 2008 
with the goal of carefully estimating the continuum level, taking 
advantage of the 1 .75 GHz bandwidth in the wide-band single- 
polarisation mode, corresponding to 4956 km s . For these ob- 
servations we used 5 antennae in the compact D-configuration. 
The final on-source integration time for this run was 5.2 hours 
(again, discarding scans with strong atmospheric phase insta- 
bilities). We detected the millimeter continuum at a level of 
S3mm=0.99±0.05 mJy at the centimeter position of the radio 
source measured by Rawlings et al. ( 1990). 



Data reduction and analysis relied on GILDAS ( Pety|2003 



The flux was calibrated against MWC349 and against reference 
quasars whose flux is monitored with the PdBI. We removed 
the continuum by assuming a point source at the position of 
3C 326 N with a flat spectrum of 1 mJy in the uv-plane. The 
data were then averaged in 71 km s~'-wide velocity channels. 

2.3. CO line emission 

We show the continuum-free CO(1-0) emission-line map of 
3C 326 N in the inset of Figure [l] integrated over 637 km s _1 
around the line centroid. The central position is slightly offset 
from the millimeter continuum position (by 0.37" x 0.80" in right 
ascension and declination, respectively), corresponding to 1/3 
to 1/6 of the beam size, respectively. Given the relatively low 
signal-to-noise ratio, this offset is not significant. 

We measured the CO(1-0) line flux in different apertures and 
in the uv plane to investigate whether the line emission may be 
spatially extended, finding a larger flux for elliptical Gaussian 
models than for a point source. This may suggest that the source 
is marginally spatially extended with a size comparable to our 
2.5"x2.1" beam. Results for the spectrum integrated over dif- 
ferent apertures are summarized in Tables[T]and[2] 

In Figure [2] we show our CO(1-0) spectrum of 3C 326 N 
integrated over a 5" aperture. The blue line shows the Gaussian 
line fit which gives a FWHM=350±100 km s _1 and an inte- 
grated emission-line flux of Ico=l-0±0.2 Jy km s _1 (over a 
5" aperture, corresponding to the slit width of the IRS spec- 
trum). The observed frequency corresponds to a redshift of 
z=0.0901±0.0001 and is similar to the observed redshift of the 
optical absorption lines within the uncertainties. 



2.4. Optical spectroscopy 

We complement these new millimeter observations of the cold 
molecular gas with publicly available optical spectroscopy from 
the Sloan Digital Sky Survey ( |York et al.|[2000| ). This spec- 
trum traces the warm ionized gas and the stellar continuum. It 
has signal-to-noise ratios of ^20-30 in the continuum (mea- 
sured at wavelengths around 5000-6000 A). 3C 326 N is sig- 
nificantly larger than the 3" diameter of the SDSS fibers, hence 
these are near-nuclear spectra including light from the central 
5 kpc. This is well matched to the 2.5" x 2.1" beam size of our 



IRAM millimeter observations, and to the 5" slit width of the 
Spitzer-IRS spectra. The data cover the spectral range between 
[OII]AA3726,3729 and [SII]AA6716,6731 emission lines and 
the Ca H+K, Mg b, Na D absorption lines. All lines are rela- 
tively broad and spectrally well resolved at the spectral resolving 
power of the SDSS of R~1800 ( [York e t al. 2000). All measured 
line properties are summarized in Table [3] While the low reso- 
lution makes it difficult to resolve certain line doublets, it does 
allow us to measure velocity dispersions and to analyze the line 
profiles. 



We used the publicly available Starlight package (Cid 
|Fernandes et al.|2005| > to fit the stellar continuum of 3C 326 N, 
and to correct our emission-line measurements for the underly- 
ing continuum emission and stellar absorption lines. Starlight 
allows to fit the optical stellar continuum of a galaxy with a lin- 
ear combination of simple stellar populations (SSPs) of different 
ages and metallicities. We used SSPs from the stellar popula- 
tion synthesis models of Br uzual & Chariot ( 2003 1 for a wide 
range of ages (from ~ 10 yrs to several 10 yrs) and metallic- 
ities (Z—0. 0002-0. 05). Wavelengths near strong emission lines 
are excluded from the fit. Formally, we give our best-fit param- 
eters. We find that 97% of the population formed 13 Gyr ago, 
with the remaining 3% formed 1.3 Gyr ago. However, visual in- 
spection of a fit with or without the intermediate-age population 
reveals no difference. Only populations with solar or superso- 
lar metallicity contribute. Adopting a Calzetti extinction law we 
find Ay =0.17 mag. In light of the unavoidable degeneracies 
and systematic uncertainties, we caution that these fits may be 
somewhat less constraining than the numerical results suggest at 
face value. This does not change the basic result, that overall, 
these fits do imply an old stellar population with age >10 Gyr, 
solar metallicity, and low extinction. We show the spectrum of 
3C 326 N along with the best fit for the most age-sensitive part 
around the 4000A spectral break in Fig. [7] Unless stated oth- 
erwise, we will in the following only refer to the continuum- 
subtracted emission lines. 

2.5. Emission-line kinematics 

After subtracting the underlying continuum emission, we fit- 
ted the emission lines with Gaussian profiles, where line 
widths, fluxes, and redshifts are free parameters, except for the 
[OII]AA3726,3729; [OIII]AA4959,5007, [NII]AA6548,6583, 
and [SH]AA6716,6731 doublets, where we required the redshifts 
and line widths in the two lines of each doublet to be identical. 
For [OIII]AA4959,5007, [NII]AA6548,6583 we also required a 
flux ratio of 1/3 between the fainter and the brighter component. 
Results for individual lines are given in Table [3] 

All lines have relatively large widths with FWHM^OO km 
s~ 1 . Redshifts and line widths are similar within a scatter of ~ 30 
km s for all lines. Redshifts are consistent with the systemic 
redshift of z=0.0900±0.0001 which we obtained from the stellar 
absorption lines. 

Careful inspection of the profiles of the relatively lu- 
minous, and spectrally well isolated [OIII]AA4959,5007 and 
[OI]AA6300,6363 lines reveals the presence of very broad com- 
ponents with widths of FWHM=1635 km s~ ! (Fig. Broad 
components of forbidden lines are not associated with typical 
nuclear broad-line regions ( Sulen tic~et al.|20 00 ) but do suggest 
that some of the ISM in 3C 326 N is kinematically strongly dis- 
turbed. The redshifts and line width of the very broad compo- 
nents are overall consistent with the wind component of the Na 
D line discussed in f|5] Equally large line widths of ionized gas 
have also been found by |Holt et al.| ( |2008] l in nearby radio galax- 
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ies showing the sig ns of AGN feedback and by Nesvadba et al.| 
(2006 2007 2008 ) in the kpc-scaled jet-driven outflows of ion- 
ized gas in high-redshift radio galaxies. 

To compensate for the difficulty of fitting broad compo- 
nents at relatively low signal-to-noise ratios with more nar- 
row, superimposed components, we required that all four lines 
have the same redshift and line widths, which yields reasonably 
good residuals. For the broad lines, we do not find a significant 
blueshift or redshift relative to the narrow components. 

Due to the large line widths, several lines are blended, 
which makes a detailed fit more difficult (namely, these 
are [OII]A3726, 3729, Ha and [NII]AA6548,6583, and 
[SII]AA67 16,6731). We did not attempt to fit multiple compo- 
nents to these lines. 

2.6. Characteristics of the Na D lines 

Figure [5 shows the strong and broad Na D absorption feature 
that we detect in our SDSS spectrum of 3C 326 N. We used the 
stellar population synthesis models described in ( 2.4 to estimate 



and remove the stellar component of the Na D line which can be 
significant in stellar populations. We also show Mg b in Figure[5] 
which is purely stellar and illustrates the accuracy of the popu- 
lation synthesis fit. 

Figure [5] shows that the Na D lines in 3C 326 N are very 
broad, so that we cannot directly measure the equivalent width 
of each component of the Na D doublet. We therefore modeled 
the lines using the atomic data compiled in [Morton ( 1991| l, al- 
lowing for changes in the turbulent velocity, velocity offsets, and 
covering fraction. Since we do not resolve individually the lines 
of the doublet, covering fraction (C f) and optical depth are de- 
generate unless the line shapes are accurately known. The lines 
are heavily blended and of insufficient S/N to have confidence in 
the line shapes. 

We find that the Na D lines are heavily saturated so that we 
can only give lower limits on the Na I column densities. Values in 
the range 4—10 x 10 14 cm~ 2 (for C/=l and 0.5 respectively) are 
most likely. The Doppler parameter of the doublet is ^800 km 



-1 



s with a most probable velocity offset of about -350 km s 
relative to the systemic redshift. For the highest column densi- 
ties, corresponding to a covering fraction C/ <0.5, the line pro- 
file begins to show a "flat core", which is not favored by the data. 
Therefore, it appears that the best fitting models favor a rela- 
tively high covering fraction. We note that the best fit Doppler 
parameter and velocity offset are not very sensitive to the range 
of columns and covering fractions explored. We will further ex- 
plore the wind properties in ( 5.1 



3. The origin of the line emission 

The line emission of interstellar gas can be powered by a number 
of different astrophysical processes like energetic photons from 
starbursts and AGN, cosmic rays, or mechanical energy injected 
into the ISM by stars, AGN, or various dynamical processes. We 
will now analyze the emission lines of the molecular and ionized 
gas in 3C 326 N to identify the physical mechanism that is giving 
rise to the luminous line emission in this galaxy. 



3. 1. Molecular emission-line diagnostics 

Emission-line ratios in the optical are a popular tool to inves- 
tigate the heating mechanism of the ionized gas (e.g., Baldwin 
[etaTl[T98T] |Veilleux & Osterbrock|[T987l ), and to differentiate 



between star-forming galaxies and AGN. Somewhat in analogy 
we will now construct a molecular diagnostic diagram based on 
the CO(1-0), pure-rotational H2 and PAH emission shown in 
Figure [6] We start with a short discussion of the astrophysical 
origin of each tracer before presenting the resulting diagram. 

Critical densities for the S(0) to S(3 ) lines are moderate (10 
to 10 4 cirT 3 in molecular gas at 500 K; Le Bourlot et al.|l999) , 
so that the lowest rotational states of H2 are populated by col- 
lisions rather than fluorescence. The first rotational lines of H2 
are therefore dominant cooling lines for molecular gas over a 
wide temperature range of ~ 10 2 ~ 3 K. They measure the power 
dissipated by all heating processes of the warm molecular gas 
(T > 100 K). This includes UV or X-ray photons produced by 
young stars in photon or X-ray dominated regions, or the dis- 
sipation of mechanical energy as proposed by |Guillard et al" 



(2009) for Stephan's Quintet. |Fe"rland et al.| ( f2008) proposed that 



cosmic rays could be the dominant energy source powering the 
extended warm intracluster H2 line in the Perseus galaxy clus- 
ter. The CO(1-0) emission-line luminosity is associated with the 
cold molecular phase, to the extent where it is commonly used to 
estimate the cold molecular gas mass (e.g., [Solomon et al.|1997) . 
Thus, if the gas traced by CO and Hi is physically associated, 
then the H2 to CO(1-0) ratio is a measure of the total heating of 
the molecular gas per unit molecular gas mass. 

PAH emission in star-forming regions is found along the 
surfaces of molecular clouds heated by UV photons, but not 



within HII regions, where PAHs are destroyed (e.g. Cesarsky 
et al.|1996 Tacconi-Garman et al.|2005 1. In addition, on galac- 



tic scales, some of the PAH emission may originate from the 
diffuse interstellar medium rather than star-forming regions, as 
argued by Draine et al. ([2007 ) based on the modeling of the dust 



emission of SINGS galaxies. In either case the PAH emission is 
powered by UV photons. Thus, we can use the ratio of PAH to 
H2 emission as an empirical measure for the contribution ofUV 
photons to the total H% heating. This is supported by the tight 
correlation between the 7.7/zm PAH feature and H2 luminosity 
found in star-forming regions (Rigopo ulou et al.|2002 Roussel 
|et al.|2007D . 

In Figure [6 we combine both line ratios into one diagnostic 
diagram. We also include the position of 3C 326 N, Stephan's 
Quintet (SQ), and of several H2 -luminous radio galaxies from 
the sample of |Ogle et al.| ( |2.010] > which have published CO mea- 
surements. We exclude central galaxies in cooling-flow clusters 
where CO line emission is extended over sizes much greater than 
the IRS slit. We also show star-forming galaxies and AGN from 
the SINGS survey, taken from Roussel et al. (2007 1. To convert 
8/im flux densities into PAH 7.7/im luminosities, we used the 
prescription of Roussel et al. (2007). We used the starburst tem- 
plate of Bran dlet al.| ( |2006] l to bring the integrated PAH fluxes 



2010) on the same 



of the radio galaxies measured by Ogle et al. 
scale. In Figure [6] 3C 326 N and the SQ shock are the most 
extreme representants of H2 luminous galaxies with the highest 
ratio of H2-to-PAH emission and H2-to-CO emission. 

3.2. What powers the Ho emission? 
3.2.1. Energetic photons 

We will now use our molecular emission-line diagnostics to 
show that UV photons produced in PDRs cannot make a major 
contribution to heating the molecular gas in 3C 326 N. 

In the diagnostic diagram shown in Figure [6] black lines 
mark the line ratios derived for PDR models to illustrate in which 
parts of the diagram the H2 heating is dominated by UV photons. 
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The H2 line emission in PDRs is derived from the calculations of 
Kauf man et al.| ( |2006| l. The total extinction of the PDR models 
is Ay =10 mag, large enough to ensure that most of the incident 
UV light is being absorbed. For the PAH emission, we assumed 
that v/ v ,8 = 0.16 X Fir, where Fir is the bolometric infrared flux, 
and 7 V ,8 the flux density at 8/zm. This conversion factor corre- 
sponds to the models of |Draine & Li|(2007| ) for a PAH-to-dust 
mass fraction qpAH — 3.55% Draine & Li (2007i, the median 



value for galaxies with Galactic metallicity in the SINGS sam- 
ple ^Draine et al.|2007) . 

The star-forming galaxies from SINGS are our control sam- 
ple with 'ordinary' line emission. They fall into the portion of 
the diagram spanned by the PDR models, as expected if most 
of their molecular line emission is powered by star formation. 
Their positions in the diagram suggest a ratio between the inten- 
sity of the UV field and the gas density of about 0.1 to 1 cm~ 3 , 
and a UV intensity several times higher than that in the solar 
neighborhood. 

In galaxies which do not fall into the portion of the diagram 
spanned by the PDR models, the molecular gas is heated by an- 
other mechanism than photons produced in star- forming regions. 
This includes 3C 326 N, 3C 424, and the intergalactic shock in 
Stephan's Quintet as the most extreme cases, as well as several 
other H2 luminous radio galaxies. Similarly, about half of the 
SINGS AGN fall into this region, albeit not at its extreme end. 

For the sake of completeness, we should note that X-ray 
heating seems also unlikely to explain the H2 line emission in the 
SINGS AGN ( |Roussel et al.p007) . Similarly, [Ogie et al.|p0T0) 
find that the observed X-ray luminosities of their sample of H2 
luminous radio galaxies are not sufficient to power the H2 line 
emission. 



3.2.2. Cosmic rays 



Ferland et al. (2008} recently proposed that heating by cosmic 
rays may explain the extended, filamentary emission of warm 
molecular hydrogen in massive galaxy clusters experiencing 
episodes of AGN feedback from their massive central galaxy. 
We discuss this excitation mechanism based on an approach that 
is broadly similar to that of Ogle et al. ( 2010) > 

The critical density of the S(0) and S(l) lines is sufficiently 
low ( Le Bourlot et al.||1999| l to allow us to safely assume that 
the J=2 and 3 levels of H2 are thermalized by collisions. The 
S(0) and S(l) line fluxes in |Ogle et al.| ( |20l"0] > correspond to an 
emission-line luminosity of 1.9 x 10 ergs . Dividing this lu - 



minosity by the H2 mass of 9.5 x 1O 8 M from |Ogle et al. ( 2010 1 
we can estimate a cooling rate through the H2 S(0) and S(l) lines 
of 3.4 x 10~ 25 ergs" 1 U^ 1 . Including the S(2) and S(3) lines in- 
creases this luminosity to 9 x 10~ 25 ergs -1 HL7 . What ionization 
rate from cosmic rays would be required to provide the corre- 
sponding amount of heat? 

Gas heating associated with the ionization of H2 by cos- 
mic rays is ~ 12 eV per ionization, including th e contribution 



from Hj" recombination and H2 re-formation (e.g. Le Petit et al. 
2006). To balance line cooling, the ionization rate per H2, L,h2, 



must be few xlO s . For such a rate, cosmic rays are the 
main destruction path of H2 molecules. The molecular gas frac- 
tion depends on the rate of ionization and gas density £7/2 /«H2, 
and the gas is molecular for n# > 10 4 cm~ 3 (see Figure 2 of 
|Ferland et al.||2008| ). At such densities, the H2 rotational states 
are thermalized up to / > 5, and the higher-J lines up to S(3) are 
also cooling lines. Dividing the H2 S(0)-S(3) luminosity by the 
energy released per ionization, we estimate the required ioniza- 



tion rate, £#2 ~5x 10~ 14 s _1 . This value is larger than what is 
inferred from observations in the Milky Way, by a factor 100 
for the diffuse interstellar medium ( Indriolo et al.|2007) , and by 
more than a factor 10 for the molecular gas within 200 pc from 
the Galactic center ( Goto et al. 2008 1. It is therefore unlikely that 



the warm molecular gas in 3C 326 N is powered by cosmic rays. 

3.2.3. Mechanical heating through shocks 

Having ruled out UV and X-ray photons as a possible mecha- 
nism to power the line emission in 3C 326 N, and having shown 
that cosmic rays are unlikely, our best remaining candidate is 
mechanical heating through shocks. This is in close analogy to 
the model of Guill ard et al.| (|2009) who studied the intergalactic 
shock in Stephan's Quintet, and found that the H2 line emis- 
sion is powered by the dissipation of kinetic energy injected by 
the interaction of two galaxies. 3C 326 N and Stephan's Quintet 
fall into similar regions of our diagnostic diagram (Figure |6}, 
suggesting their line emission is powered by a similar physi- 
cal mechanism. In §Sj4]and[6]we will demonstrate that the radio 
source is the only viable energy source which can explain all of 
the phenomenology observed in 3C 326 N, and present a phys- 
ical framework through which the energy injection by the radio 
source and the molecular line emission can be related. 



3.3. What powers the optical line emission': 



In the previous section we presented new, 'molecular' diagnos- 
tics to argue that the luminous H2 emission in 3C 326 N (as well 
as other H2 radio galaxies) may be powered by the dissipation of 
mechanical energy in the interstellar medium of the host galaxy. 
It may be illustrative to use our SDSS spectra to compare this 
result with the optical diagnostics. 

Table|3]lists the fluxes measured for various optical emission 
lines in 3C 326 N. Comparison with the classical BPT-diagrams 
([Baldwin et al.|1981||Veilleux & Osterbrock|1987||Kewley et al 



2006) shows that 3C 326 N falls within the LINER ( [Heckman 
1980) part of the diagrams, similar to most H2 luminous radio 
galaxies of O gle et al.|(|2010|l and also the galaxy-wide shock in 
Stephan's Quintet ( |Xu et al.|2003) . 

The excitation mechanism of the atomic lines in LINERs is 
a long-standing issue in the literature where some studies favor 



photo-ionization by the AGN (e.g., Veilleux & Osterbrock 1987[l 



while other studies propose shock excitation (e.g., |Clark et al. 



1998 Dopita et al. 


1997 Monreal-Ibero et al.|2006|l.|Dopita & 


Sutherland ( 1995 


1996) show that the spectral characteristics 



of LINERS can be modeled with excitation from fast, radiative 
shocks without the emission from the radiative precursor of the 
shock. The radiative precursor of a fast shock moving into low- 
density gas adds an emission component which has the spectral 
characteristics of photoionized gas with a high ionization param- 
eter. However, the line ratios in 3C 326 N, in particular the ratio 
between the [NeII]/[NeIII] lines in the mid-infrared, suggest a 
low ionization parameter. This apparent contradiction can be re- 
solved if a clumpy or filamentary gas distribution enhances the 
'leaking' of hard photons along sight lines where the photons 
can escape without interacting with the pre-shock gas. As in the 
shock of Stephan's Quintet, most of the gas in 3C 326 N is in H2 
at densities nu > 10 3 cm -3 much higher than the density of the 
ionized gas (Guillard et al. 2009), in good agreement with this 



Quantitatively, using the models of Allen et al. ( 2008[ ) for 
3C 326 N, we find that the ratios between the main optical lines 
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in Tabled are close to those of a 250 km s 1 shock moving into 
ionized gas with a 1 cm~ 3 preshock density and a low magnetic 
parameter, B/y/n <1 /iG cm -3 / 2 . Observed and modelled line 
ratios are compared in Table [4] This is not a unique solution, in 
particular, the line ratios are also compatible with lower densities 
for similar shock velocities. A detailed analysis of the alternative 
models (and model parameters, as done in |Dopita et al.f l997 ) 
for the interaction between radio jet and interstellar medium in 
M87 is beyond the scope of this paper. 

We will merely use these models to determine the bolometric 
luminosity of the ionized gas from the observed H/3 luminosity. 
Using the output of |Allen et al. (2008 ) for shock velocities below 
about 300 km s _1 we can calculate the ratio of H/3 luminosity, 
-2//^, to the total luminosity produced by the shock. We find 
that J£ s ~ 200 ± 50J^// p . This bolometric correction is slightly 
larger for magnetic parameters below the equipartition value. 
This large correction factor is due to the much larger brightness 
of the UV emission lines relative to the optical line emission, 
which we did not observe. With the measured H/3 flux and using 
this model at face value, we find Jz? s ~ 1 x 10 43 ergs s _1 , a factor 
~10 more than the total H2 line luminosity in the S(0)-S(7) ro- 
tational lines ( Ogle et al.|2010 l. This corresponds to the isobaric 
cooling of gas at temperatures of ~ 10 6 K for a mass flow of 
order of a few 100 M yr _1 . This mass flow is very large, and 



Ogle et al.| (|2010| find an X-ray luminosity of 4 x 10 40 erg 



we will discuss in S6.3 how this result may be interpreted as a 
repeated mass cycling between gas phases driven by interactions 
between the radio jet and the multiphase interstellar medium of 
the host. Obviously, due to the uncertainties in the measurement 
and modeling all of these values have uncertainties of factors of 
a few. 

It is well possible that not all of the optical line emission 
is excited mechanically. The dynamical interaction between gas 
phases must include turbulent mixing between the cold and hot 
gas, which produces extreme UV radiation along the surfaces of 
cold clouds. This radiation is an additional, and maybe signifi- 
cant, source of ionization powered by the thermal energy of the 
cocoon. ICrawford & Fabian (1992 ) show that this mechanism 



could account for parts of the line emission in ionized gas in 
cooling flows, consistent with observed line ratios. In this case, 
the bolometric correction with respect to H/3 depends on the ef- 
fective temperature of the gas after mixing and could be sig- 
nificantly smaller. With regard of this uncertainty on the radia- 
tive excitation, we will in the following adopt a fiducial value 
of 1 x 10 43 erg s . Again, the uncertainties of this estimate are 
likely factors of a few. 



3.4. Photoionization by stars and the AGN 

In the above estimates we assumed that all of the optical line 
emission is due to shocks, and that other processes like pho- 
toionization from star formation or the AGN can be neglected. 
With the low star-formation rate and faint AGN X-ray emission 
from 3C 326 N this galaxy is ideally suited to set stringent limits 
on the contribution of these processes to the observed line emis- 
sion. For example, the observed star-formation rate of < 0.07 
M yr~ ! would imply an Ha luminosity „£?(Ha)~ 1 x 10 40 erg 
s _1 (corresponding to a H/3 luminosity of few x 10 39 erg s _1 ) 



s 1 for 3C 326 N, whereas our [OIH]A5007 measurements in- 
dicate an [OIII]A5007 luminosity of log«5f [0///] = 7 x 10 40 erg 



s (including the extinction correction of a factor 2). Heckman 



et al. ( 2004) > find that the bolometric luminosity of quasars 
scales with [OIII]A5007 emission-line luminosity as «Sff, ; ~ 
3500^j o///]A5007 . For 3C 326 N the observed [OIII]A5007 flux 
implies a bolometric luminosity of few x 10 44 erg s _1 , three or- 
ders of magnitudes greater than that estimated from the X-ray 
measurement. This indicates that the AGN photoionization is 
unlikely to be the dominant mechanism in exciting the optical 
line emission in 3C326 N. 



4. Mass and energy budgets 

The relative mass budgets of warm and cold molecular gas 
and the warm ionized gas provide important constraints on the 
physical and astrophysical conditions of a galaxy. By 'physi- 
cal conditions' we refer to the gas properties, which relate di- 
rectly to the observations. By 'astrophysical conditions' we re- 
fer to the galaxy properties and the mechanism which is causing 
the observed gas properties. Typically, in gas-rich, actively star- 
forming galaxies the amount of cold molecular gas exceeds the 
masses of warm molecular and ionized gas by factors 10—100 



(see Roussel et al. 2007 1 Higdon et al. 



2006 for samples of 



nearby galaxies and ULIRGs, respectively). We will now show 
that this ratio is much smaller for 3C326 N. All masses, lumi- 
nosities, and kinetic energies are summarized in Table [5] 



4. 1. Molecular gas mass 

4.1 .1 . Direct estimate of the luminosity and mass of warm 
molecular gas 

The mass of warm molecular gas in 3C 326 has been estimated 
by |Ogle et aL] p007] |20T0| > by fitting the H 2 S(0) to S(7) rota- 
tional line fluxes with 2 or 3 components at different tempera- 
tures where H2 excitation and ortho-to-para ratios are assumed 
to be thermalized. This yields a mass of 1.0 x 10 9 M . 

In the present analysis, we associate the H2 emission with the 
dissipation of kinetic energy in the molecular gas. Therefore we 
follow a different approach to estimate the mass and luminosity 
of warm molecular gas. As in Guillar cTet al.| l |2009) , we model 
the dissipation process with magnetic shocks in molecular gas, 
which maximizes the H2 luminosity per total emitted power. In 
this sense, our results represent a lower limit to the dissipated 
energy required to account for the observed H2 luminosity. 

Each shock model includes a range of gas temperatures 
which depend on the shock velocity, pre-shock gas density and 
intensity of the magnetic field . We use the grid of models pre- 
sented by Guillard et al. (2009 1 for proton densities of n# = 10 3 



cm -3 and 10 4 * cm~ J , respectively, an initial ortho-to-para ratio 
of 3 and a magnetic parameter, B / y/riH = 1 /xG cm~ 3 / 2 . 

The shock velocity is the only parameter that we allow to 
vary. A combination of three shocks is required to match the 
emission in all 8 H2 lines, S(0) to S(7). These fits provide a scal- 
ing factor for each of the three shocks which represents a mass 
for a continuous star-form ation history and using the models of fl w, the am ount of gas traversing the shocks per unit time. As 
( |2003| ._This corresponds to ~ 5% of the mea- discussed in |Guillard etaL] (|2009|), this fit is not unique but the 



Bruzual & Chariot 



sured Ha luminosity of 3C 326 N. 

We can also rule out a dominant role of the AGN in pho- 
toionizing the gas, for the simple reason that the [OIII]A5007 
luminosity alone exceeds the luminosity emitted in the X-ray. 



proposed solution may be used to quantify the relevant range of 
shock velocities, and to estimate the warm gas masses by multi- 
plying the mass flows with the gas cooling time, and the shock 
luminosities by integrating over all cooling lines. 
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Excitation diagrams for models with two different gas densi- 
ties, 10 3 cm~ 3 and 10 4 cm~ 3 , respectively, are shown in Fig. [3] 
The corresponding H2 line fluxes are listed in Table [6] They are 
smaller than the kinetic energy fluxes one may compute from 
the mass flow and shock velocity because some of the energy is 
transferred to the magnetic field. The luminosities of the mid- 
infrared H2 lines are close to the bolometric values obtained by 
integrating the emission in all lines because the H2 rotational 
lines are the main cooling lines of magnetic shocks. Summing 
over the three shocks we estimate a total luminosity of 10 42 erg 
s _1 for the molecular Hydrogen. 

Table [7] lists the gas cooling times and H2 masses down to 
temperatures of T=150 K for each of these models. The total 
mass we obtain for gas at temperatures larger than T=1 50 K, 
1.3 — 2.7 x 10 9 M Q , is slightly larger than that derived by Ogle 



et al. (2010 1, and we will in the following assume a fiducial 
molecular gas mass of 2 x 10 9 M . The small difference be- 



tween the estimates of |Ogle et al. (2010 1 and our results is due 
to different ortho-to-para ratios and moreover, in our lowest- 
density models, the S(0) and S(l) lines are not fully thermalized. 



4.1 .2. CO(1-0) luminosity and total molecular gas mass 

The CO emission-line luminosity is often used as an empirical 
measure of the mass of cold molecular gas. In { 2.3 we estimated 
an integrated CO(l-O) emission-line flux of Icq=1.0±0.2 Jy 
km s _1 from our millimeter spectroscopy at the IRAM Plateau 
de Bure Interferometer, extracted from a 5" aperture. Using 
Equation 3 of Solo mon et al.| ( |T997] > we translate this value into 
a CO(l— 0) emission-line luminosity of Leo' — 3.8x 10 8 K km 
s _1 pc 2 at a redshift of z=0.090. Applying a Galactic H2-to- 
CO conversion factor of 4.6M /(K km s~ ! pc 2 ) as determined 
for gas in the molecular ring of the Milky Way ( Solomon et al. 
1992), this would correspond to a mass of cold molecular gas of 
M^° d N = 1 .5 x 10 9 M . We will further discuss the appropriate- 
ness of this H2 conversion and related uncertainties in f|7] 

Comparison with [4.1.1 shows that the amount of warm 
molecular gas measured directly from the mid-infrared lines is 
similar to the H2 molecular mass inferred from the CO(1-0) 
flux estimated from our PdBI observations. This is a highly un- 
usual finding compared to star-forming galaxies where the ratio 
of warm to cold molecular gas mass is of order 10~ ! to 10~ 2 
( |Roussel et al.|2007[|HTgdon et al.|2006| for the same CO-to-H 2 
conversion factor. However, it is not very different from the ra- 
tio of ~ 0.3 found in molecular clouds near the Galactic center 
( Rodriguez-Fernand ez et al.|2001| l. 

The Galactic center may be an adequate nearby analog for 
the properties of the molecular gas in 3C 326 N. In the Galactic 
center cold dust temperatures indicate that the gas cannot be 
heated by UV p hotons ( Lis et al.|2001 1, but possibly by shocks 
or cosmic rays (jRodrfguez -Fernandez et al. 2001 ; Yusef-Zadeh 



et al. 



Goto et al. 



(2008 



2007 1. The Ht observations presented by 
imply an H2 ionization rate lower than the value required to ac- 
count for the temperature of the warm molecular gas of ~ 150 K 
(Rodriguez-Fernandez et al. 2001) with heating by cosmic rays, 



C,m ~ lQ-^s- 1 |Yusef-Zadeh et al.|2007 1. 

From detailed studies of several molecular species in the 
Galactic center, |Lis et al.| ( |2001| l; |Huettemeister et al.| (1998 ) 
conclude that all of the molecular gas in this environment may 
be warm at temperatures above T=50 K. By analogy we cau- 
tion that the usual distinction between cold and warm gas mea- 
sured through CO line emission and infrared H2 lines, respec- 
tively may not apply to 3C 326 N. Most of the CO line emis- 
sion in 3C 326 N could in fact be associated with the warm gas. 



Observations of higher-J transitions are necessary to test this hy- 
pothesis. We used the RADEX LVG code (van der Tak et al. 
2007 ) to verify that the intensity of the CO( 1 -0) line is relatively 
insentitive to the gas temperature for a given CO column density. 
For example, in a temperature range of 10-100 K, the intensity 
changes by about a factor 2 (see also Lada & Fich|1996 i. 

Obviously, given these considerations, the standard conver- 
sion factor may not apply at all. With these caveats in mind, we 
will in the following assume that the two tracers do not represent 
the same gas, which would correspond to a total gas mass of at 
most 3.5 x 10 9 M n . 



4.1.3. Warm ionized gas mass 

We can also estimate the mass of ionized gas in 3C 326 N from 
our SDSS spectra. Assuming case B recombination and follow- 
ing |Osterbrock| (jT989j we can estimate an ionized gas mass, 
Mhh from the Ha emission-line luminosity, 5£na, by setting 



Mhii = 



h V Ha &Ha 



= 9.73 x 10° jSf ffo ,4i n;l m M Q , 



(1) 



where m p is the proton mass, h is Planck's constant, a e J^ and 
VHa are me effective recombination cross section and the fre- 
quency of Ha, respectively. n e ,ioo is the electron density in units 
of 100 cm~ 3 . The Ha luminosity, J?// a ,4i, is given in units of 
10 41 ergs" 1 . 

In order to estimate extinction-corrected, intrinsic Ha lu- 
minosities, we measure the Ha and H/3 fluxes from our 
SDSS spectra and compare with the expected Balmer decre- 
ment of Ha/H)3=2.9. We find A™ N =0.% mag. This implies 

an extinction-corrected Ha luminosity of Jz? (//a) 32 , 6 N — 2.8 x 
10 41 ergs" 1 . 

We also need to constrain the electron densities in the op- 
tical emission-line gas, which we estimate from the line ratios 
of the density-sensitive [SII]AA67 16,6731 emission line dou- 
blet. Fitting each line of the doublet with Gaussian profiles, we 
find a line ratio of R 32&v = 1.3. This is near the low-density 
limit, and assuming an electron temperature of 10 4 K, these 
line ratios correspond to electron densities n 326 =130 cm -3 . With 
these extinction-corrected Ha luminosities and electron densi- 
ties, and using Equation [T[ we estimate ionized gas masses of 



M 



326 
Hll 



2 < 10 M . 



4.2. Kinetic energy provided by the radio source 

3C 326 has a powerful, large radio source of Mpc size and re- 
sides in a relatively low-density environment compared to pow- 
erful radio galaxies generally (Stocke et al. 1979 Willis & Strom 
1978[). Estimating the intrinsic properties of radio jets, in partic- 



ular their kinetic power and lifetimes, remains a challenge. The 
synchrotron emissivity of the radio jet depends on the intrinsic 
jet power and the ambient conditions, and the jet kinetic power 
is therefore not easily derived. Estimates based on theoretical ar- 
guments suggest factors of ~10— 1000 between radio luminosity 
and kinetic power (e.g. |De Young|1993 Bicknell et al.|1997) . 

Given these uncertainties, it may be best to use empirical 
estimates of the jet kinetic power. Birzan et al. ([2004 2008} es- 
timated the kinetic power of radio sources by comparing the nec- 
essary energy to inflate X-ray cavities in galaxy clusters with the 
luminosities of the radio sources that are inflating them. They 
find a factor ^100 between kinetic and bolometric radio lumi- 
nosity, with a large upward scatter of up to factors of a few 1000, 



8 



Nesvadba et aL: The multiphase gas in a H2 luminous radio galaxy 



perhaps suggesting that the correlation is in fact a lower enve- 
lope. Using their scaling with the radio luminosity measured at 
327 MHz, and measured radio flux, we find kinetic luminosities 
4 x 10 44 erg s~ 1 . For fluxes measured at 1 .4 GHz, 

Using the estimate of 



Of ^,327 



we find 



1.3 x 10 45 erg s _1 . 



Merloni & Heinz (2007) instead, which relies on the 5-GHz core 



10 44 erg s- 1 



Willis & Strom 



(1978} 



radio power, we find 2 x 

suggest a total energy content of 6 x 10 erg in the radio lobes, 
which would correspond to a kinetic luminosity of fewx 10 44-45 
erg s for a fiducial lifetime of the radio source of 10 7 ~ 8 yrs, 
in good agreement with the previous estimates. The estimates 
derived with each method change, reflecting the uncertainties 
of each approach, but the overall result of this section holds, 
namely, that the kinetic luminosity is of order fewx 10 44-4 -" 1 



erg 



A closely related quantity is the lifetime of the radio source. 
Jet lifetimes may either be estimated based on radio spectral in- 
dices or estimates of the velocity with which the jet expands, 
the former typically giving significantly smaller values. For 
3C 326 specifically, Willis & Strom ( 1978] > carefully investigated 
the multi-wavelength radio properties finding spectral ages of 
1 — 2 x 10 yrs. Using their estimate of the Alfven speed and 
size estimate with our cosmology would yield a kinetic age of 
6 x 10 7 yrs for most parts of the radio source, although a region 
that is somewhat spatially offset and has a different polarization 
angle, may be as old as ~ 2 x 10 yrs ( |Willis & Stroml|l978| . 
It is unclear whether this is due to an extended activity period, 
several activity outbreaks, or the second galaxy of the system, 
3C 326 S, which also has a radio core. At any rate, the detection 
of millimeter continuum emission from 3C 326 N and 3C 326 S 
({2.2 see also Fig. [I} implies on-going activity in both nuclei. 



4.3. Kinetic energy of the gas 

We will now give rough estimates of the related kinetic energies 
in each gas phase i. The total kinetic energy estimated from the 
line widths is given by E r . to t = 3/2 £ m,- of. m; is the gas mass in 
each phase, and (J,- is the velocity dispersion for the ensemble of 
clouds derived from the widths of the emission lines (and given 
in Tables [T] and [3] for CO and the optical lines, respectively), 
where we set a = FWHM/2.355. 

From the optical and millimeter spectroscopy we have di- 
rect measurements of the line width, which is FWHM^ 600 km 



s _1 for the optical lines, and FWHM— 350 km s _1 for the CO 
lines. The rotational H2 lines are not spectrally resolved with 
IRS, yielding an upper limit of < FWHM~ 2500 km s _1 , which 
is not very constraining. Near-infrared observations of the ro- 
vibrational lines would be an important test to infer the dynami- 
cal state of the warm molecular gas. The line widths of CO and 
the optical lines, which trace gas that is colder, and gas that is 
warmer than the warm H2, respectively, are not very different 
(Figure [2}. We therefore expect that the line widths of warm H2 
will be in the same range. 

Since these spectra are integrated, the line widths will be af- 
fected by large-scale motion within the potential of the galaxy. 
|Holt et aL] ( |2008) > obtained longslit spectra of optical emission 
lines in 16 nearby powerful radio galaxies, finding that the kine- 
matics are often very irregular and do not appear dominated by 
rotation. Where they found regular velocity gradients resembling 
rotation curves, their data suggest that within radii of 2.5 kpc 
(corresponding to our aperture) we may sample a range of rota- 
tional velocities less than 100—150 km s _1 . 



For a measured line width of FWHM^500— 600 km s a 
velocity gradient of ~ 150 km s subtracted in quadrature will 
lead to a negligible correction of ~25 km s (35 km s _1 for 
the measured CO line width of FWHM=350 km s _1 ). Therefore 
we do not believe that gravitational motion will have a large im- 
pact on our measurements. With the velocity and mass estimates 
given above and in Sj4l respectively, we obtain a total kinetic en- 
ergy of E^ = 3x lCr 7 erg for the cold and warm molecular gas. 
The ionized gas has a much smaller mass and a negligible kinetic 
energy of order 5 x 10 55 erg. Kinetic energies are summarized in 
Table |5] 



5. A wind without a starburst 

5. 1. Mass and energy loss rate of the neutral wind 

In attempting to estimate the characteristics of the multiphase 

Iwe 



medium, absorption lines can play an important role. In 1 2.6 



described the detection of a significant interstellar component to 
the Na D absorption line in 3C 326 N with a systematic velocity 
offset to the blue and a pronounced blue wing. This is very sim- 
ilar to what is frequently found in starburst galaxies (Heckman 
|et al.|2000t|Martn|2005||2006) , where blueshifted Na D absorp- 
tion is commonly interpreted as strong evidence that galaxies 
with intense star formation are driving energetic outflows. 



Similarly for radio-loud AGN, |Morganti et al. ([2005) give 
compelling evidence for outflows of neutral material based on 
studies of HI absorption line profiles at 1.4 GHz. They iden- 
tify pronounced blueshifted components with velocities of up to 
~ 2000 km s _1 in a number of galaxies, and a clear excess of 
blueshifted relative to redshifted material. They emphasize that 
they find blueshifted material in all radio galaxies with suffi- 
ciently deep HI spectroscopy, suggesting that outflows of neutral 
gas may be common amongst powerful radio-loud AGN. Three 
galaxies observed with Spitzer-IRS have also been observed in 
( 2005| l, and interestingly the two which 



HI by Morganti et al. 



have the more pronounced HI absorption have also luminous H2 
line emission (Ogl e et al.|2010 Haas et al.|2005b[ l, as expected 
if outflows and H2 line emission are physically related. The rel- 
ative velocity of the blue wing of the Na D line suggests a termi- 
nal velocity of ~ — 1800 km s~ ! for 3C 326 N, within the range 



found by Morganti et al. (2005 1. 



To estimate the energy and mass loss rates of the outflow, we 
need to constrain the associated column density of neutral gas. 
The relationship between the column density of Hydrogen, N#, 
and that of Sodium, Njvo. is determined by the abundance of Na 
relative to H and the ionization correction for Na (Nal and Nail). 
If, as is common for starburst-driven winds (e.g., Heckman et al. 



2000), we assume a solar abundance ratio, a depletion factor of 
10 ( |Morton|1975| > or perhaps more (Phillip s"et al.|1984| >, an ion- 
ization correction of a factor of 10, and a Nal column density 
of 10 14 cm~ 2 , we find a total H column density of N#~10 2L7 



cm 2 . To be conservative, we will adopt N//^10 



21 -2 

cm as a 

fiducial value. This opens the possibility that 3C 326 N may have 
a significant component of warm neutral gas, but given the ab- 
sence of HI spectroscopy, the exact amount of this gas is difficult 
to quantify, in particular, since parts of the Na I could be associ- 
ated with outflowing molecular gas. 

To directly associate the kinematics of the gas with an out- 
flow rate of mass and energy, we assume a simple model of a 
mass conserving flow with a constant velocity, which extends 
from some minimal radius to infinity. This gives mass outflow 
rates of: 
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In addition, the terminal velocity of —1800 km s is ~3- 
4 x greater than terminal velocities typically fou nd in starburst- 



-M yr-S (2) 



and energy loss rates of 



1N H 



io 42 - 5 c 



N H 



Q 



An lkpc 10 21 cm- 2 v 350kms 



-Y ergs- 1 ^) 



respectively, where Cf is the covering fraction, J| is the opening 
angle, v is the outflow velocity. Covering fraction, opening angle 
and column density are degenerate. However, relevant for our 
analysis is the product of the three quantities, so that this does 
not influence our results. But note that our analysis can only be 
accurate at an order-of-magnitude level. 

We have argued that the NaD lines are likely associated with 
a significant (but difficult to quantify without direct HI observa- 
tions) reservoir of warm atomic gas. The size of the (marginally) 
resolved CO emission (2.5"x2.1") is likely similar to the size of 
the Na D absorbing region. If we take the geometrical mean of 
the half beam width as the entrainment radius, ^1.8 kpc, allow- 
ing for a covering fraction of 0.5 to 1, and corresponding column 
densities of 10 21 " 3- 21-7 cm -2 , an opening angle of the outflow of 
n, and an outflow velocity of 350 km s _1 corresponding to the 
offset velocity found with our line fit (and ignoring the Doppler 
parameter which suggests higher velocities), we find a mass out- 
flow rate of about 30-40 Mq yr _1 and an energy loss rate of 
~10 42 - 43 ergs" 1 . 

In starburst galaxies, the energy loss rates estimated with 
the same method and assumptions correspond only to a few 



percent of the injected mechanical power (e.g. Heckman et al. 
2000 IMartin 2005). If the same factor approximately applies 



for 3C 326 N, we w ould expect a mechanical power of at least 



few xlO 43 erg s" 1 (jj4.2i 



5.2. Outflow energetics 

The outflow properties of 3C 326 N are very reminiscent of 
what is observed in local starburst galaxies (e.g. Heckman et al. 



2000 



Martin 



2005 ), which have column densities of <■ 



fewxl0 zJ 

cm - ', mass and energy outflow rates of a few to a few 10s 
M Q yr _1 and 10 43 erg s _1 , respectively, and reach maximal 
blueshifts of ^400—600 km s _1 . AGN and starburst activity of- 
ten coincide making it difficult to uniquely identify the energy 
source that is driving the wind. For 3C 326 N this is not the 
case. |Ogle et al.| ( f2007| l estimated an upper limit of SFR<0.07 
M Q yr _1 for the total star formation rate. Assuming that all 
of this star formation occurs within regions covered by the 3" 
SDSS fibre, this would correspond to a star-formation intensity 
of SFI< 0.004 Mq yr~' kpc -2 . This is more than an order of 
magnitude lower than the SFI^ 0.1 M Q yr _1 kpc~ 2 threshold 
found by Heckman ( 2003) 1 for starbursts that drive winds. We 
come to a similar conclusion when comparing with the super- 
nova rate expect ed for constant star forma tion with SFR= 0.07 



M Q yr 1 . Using Bruzual & Chariot ( 2003| > we find a supernova 



rate of 3.5 x 10~ 4 yr~\ Assuming that all of the "canonical" en- 
ergy release of a supernova, 10 51 erg, will thermalize, we find an 
energy injection rate of 10 40 erg s , about three orders of mag- 
nitude lower than the mechanical power of the wind in 3C 326 N 
or the line emission of the warm ionized and molecular gas. For 
more realistic thermalization efficiencies of few tens of percent 
( Strickland & Heckman|2009 1, the energy injection from super- 
novae will be even lower by factors of a few. 



driven winds (which are of order few 100 km s~ Heckma n et al. 
2000; Martin 2005, 2006) and similar to some of the most pow- 
erful starbursts (which have star-formation rates orders of mag- 
nitude higher than 3C 326). Both arguments indicate that the 
outflow in 3C 326 N is not related to star formation. This ve- 
locity is also significantly larger than what we may expect from 
velocities due to a possible interaction with the nearby galaxy 
3C 326 S (see the discussion in |Ogle et aL] 2007[ ). The only 



plausible candidate driving this outflow is a radio-loud AGN 

6. Energy and momentum exchange in the 
multiphase cocoon 

We will now use our observational results to construct a sce- 
nario for the interaction between jet and interstellar medium of 
the host galaxy, which explicitly includes the molecular gas. We 
will argue that the kinematics and emission-line luminosities of 
3C 326 N are most likely related to the energy injected by the 
radio jet, which is being dissipated by the multiphase interstellar 
medium of the host galaxy. We describe and quantify the energy 
flow and associated timescales. 

6. 1. Energy, mass and momentum flow 

Much theoretical effort has been dedicated to describing the 
interactions of radio jets with the ambient medium (for early 
studies see, e.g., |Scheuer 1974| Begelman & Cioffi| [l989). 
Observational evidence that jets may have a profound influence 
on the interstellar medium of their host galaxies has been known 
for at least two decades (e.g., van Breugel et al.| 1985| Pedlar 
|et al.|1985|rTadhunter|199"Tj|Eales & Rawlings|1993| l. The 'co- 
coon model' describes interactions between radio jet and ambi- 
ent gas in form of a 'waste energy basket' ( Scheuer||1974| l of 
hot, low-density, but high-pressure material that surrounds the 
thin relativistic jet. This 'cocoon' will expand into the ambient 
gas and may entrain matter ablated from denser clouds of colder 
material (e.g., Begelman & Cioffi 1989), drastically enhancing 
the efficiency with which the jet interacts with the gas of the 
galaxy compared to simple interactions along the jet working 
surface ('dentist's drill', Scheuer 1982). We will now extend this 
scenario by explicitly taking into account the energy, mass, and 
momentum exchange between different gas phases which result 
from the injection of mechanical energy by the radio source. This 
scenario represents the synthesis of our observational results dis- 
cussed in §Sj3]to|5] 

Fig. [8] illustrates our basic scenario of the energy and mo- 
mentum flow in the cocoon. Part of the mechanical energy in- 
jected by the radio jet is ultimately transformed into thermal en- 
ergy of the warm atomic and molecular gas, giving rise to the 
observed line emission. A part of the jet kinetic power is also 
translated into thermal and bulk kinetic energy of the hot co- 
coon plasma (first box in Fig. [8j. This plasma is strongly over- 
pressurized and expands on a timescale shorter than its cooling 
time, thus driving a net outflow of multiphase gas with outflow 
rates of a few xlO M Q yr -1 (Sj5j. In addition, the ablation of 
cloud material by the outflowing hot medium may replenish the 
hot medium of the cocoon, which would help maintaining a high 
pressure as the cocoon expands, and enhance its lifetime. The en- 
ergy cascade causes a momentum transfer from large-scale bulk 
motion to turbulent motion on smaller scales. Dynamical inter- 
actions between the different gas phases drive fragmentation of 



10 



Nesvadba et aL: The multiphase gas in a H2 luminous radio galaxy 



the molecular clouds and turbulent motions between and within 
individual fragments (second box in Pig . [8]> . This may lead to en- 
trainment of parts of the warm and cold medium in the hot wind 
as described for starburst-driven winds (Hec kman et al.|2000| l. 

This interaction between the molecular gas and the out- 
flowing plasma creates a physical environment not very differ- 
ent from that produced by the galaxy-wide shock in Stephan's 

|Appleton et al.|2006) . Similar to 



Quintet (Guillard et al. 



the analysis of |Guillard et a 



2009 



(2009), 



we postulate that turbu- 
lent motion in the cocoon drives shocks with velocities that de- 
pend on the gas density. In a multiphase medium, differences in 
gas densities between the cold and hot gas are up to factors of 
10 6 . For a given ram pressure, this translates into a factor 10 in 
shock velocity. The shocks maintain the amount of warm molec- 
ular and ionized gas that are necessary to explain the observed 
luminous line emission of the molecular and ionized gas (third 
box in Fig. 5J. Due to the high density, shocks driven into mag- 
netized mo ecular gas will be slow, with velocities of a few 10s 
of km s , so that H2 will not be destroyed, but becomes a main 
coolant. This is in agreement with our molecular diagnostic dia- 
gram in Fig.|6]and with the pure rotational H2 line ratios, which 
are consistent with excitation through slow shocks. 

The gas cooling times for these shocks listed in Table [7] are 
very short, of order 10 4 yrs. In order to maintain the gas at the ob- 
served temperatures, the cocoon must inject energy over similar 
timescales. This is important, since these timescales are signifi- 
cantly shorter than the free-fall times of self-gravitating molecu- 
lar clouds. This is evidence that the gas cannot form gravitation- 
ally bound structures and stars, simply because it is continually 
being stirred up by mechanical interactions in the multiphase gas 
which are ultimately powered by the radio source. We will dis- 
cuss in §6. 3 1 that the AGN may maintain such conditions over 
significant timescales (10 yrs), and discuss astrophysical im- 
plications in S]7] 

6.2. Efficiency of the energy transfer 

Simple energy conservation implies that our scenario for pow- 
ering the line emission and the outflow through mechanical in- 
teractions is only realistic if the efficiency of the power transfer 
in each step is < 1 . We hav e estimates for the jet kinetic power 
(>few xlO 44 erg s _1 , (4.2 1 and the emission-line luminosity o:" 
H2 and HII (~ 10 43 erg s _1 ). The mechanical power necessary 
to drive the outflow is similar to the emission-line luminosity at 
an order-of-magnitude level ((j5}. Taken at face value, these es- 
timates indicate that at least ~10% of the jet kinetic luminosity 
is deposited within the cocoon, which shows that our scenario 
of jet-powered line emission is indeed energetically plausible. 
The efficiency of the cocoon in driving an outflow - the ratio be- 
tween the mechanical power carried by the outflow to the power 
thermalized in the cocoon - may well be about 1 , comparable to 
that derived from observations of starburst-driven outflows (e.g. 
M82 Strickland & Heckman 2009 ). The luminous line emission 
of 3C 326 N implies that a power comparable to that of the bulk 
outflow is radiated away due to dissipation in the interstellar 
medium. 



6.3. Dissipation time and length of the H 2 -luminous phase 

If the radiation of luminous optical and infrared line emission is 
due to the dissipation of the kinetic energy of the warm and cold 
interstellar medium, then we can roughly estimate the dissipa- 
tion time of the kinetic energy of the gas, simply by relating the 



measured kinetic energy and the emission-line luminosity. We 
will in the following estimate this timescale, and compare with 
the characteristic timescales of jet activity. 

In (3.3 we found that the ionized gas, although negligible 
in the mass budget, dominates the radiative energy budget of the 
warm and cold interstellar medium with a bolometric luminosity 
of ~ 1 x 10 43 erg s _1 , whereas the molec ular gas emits a lower 
luminosity of 1 x 10 42 erg s _1 . 



In S4.3 



we estimated that the 
kinetic energy of the molecular gas is most likely ~3x 10 57 
erg, about two orders of magnitude more than that of the warm 
ionized gas (see also Table [5]), and that it is not likely that this 
energy budget will be dominated by rotation. 

We will now constrain the range of plausible dissipation 
times, assuming that the UV-optical emission of the ionized gas 
is powered by the turbulent kinetic energy, and fully participates 
in the dissipation process. In this case, 



3/2 {M ion af on +M mul al ol ) 



mol 



10 7 yr 



(4) 



But the dissipation time could be significantly larger if the 
UV-optical lines are not entirely powered by the turbulent en- 
ergy. The molecular gas could be in a region where the hot 
plasma has already been accelerated, so that the kinetic energy 
of the hot wind would power the line emission through direct in- 
teractions between the molecular gas and the surrounding warm 
and hot gas. In this case, the bulk kinetic energy of the wind 
would be powering the line emission, not the turbulence of the 
molecular gas. In the extreme case that the luminosity of the 
ionized gas is fully dominated by the kinetic power from the out- 
flow, it does not contribute to the dissipation and the timescales 
simply follow from setting 



3/2 {M mol al ol ) 



10V 



(5) 



(Note that we included He into our mass budget for both time 
estimates.) These findings have two interesting consequences. 
First, the dissipation time is much longer than the cooling time 



predicted by our shock models for the molecular gas ((4.1.1 
which is of order 10 4 yrs. This may imply that the turbulent en 
vironment of the cocoon feeds a mass and energy cycle similar 
to that described by Guillard et al. (2009) for Stephan's Quintet, 
where the gas undergoes many episodes of heating and cooling 
on "microscopic" scales, effectively maintaining an equilibrium 
between different gas phases on macroscopic scales. 

The large mass flow of a few 100 M yr _1 (which we esti- 
mated in 



3.3 



from the H^ luminosity) is a natural outcome of 
this mass cycle. It is too large to be accounted for by the large- 
scale bow shock as the cocoon expands through the galaxy, but 
may plausibly be produced by shocks that are locally generated, 
as the molecular gas fragments move relative to the low-density, 
hot medium. Shock velocities of order 250 km s inferred from 



the optical line ratios (( 3.3 1 are subsonic with respect to the hot 
gas, but supersonic relative to the warm gas (T < 10 6 K). Such a 
cycling between warm and cold gas phases may also be impor- 
tant for the exchange of momentum between the hot plasma and 
the H2 gas and thereby for the entrainment of gas into the flow. 

Second, a dissipation time of ~ 10 7 ~ 8 yrs is signific ant co m- 
pared to the lifetime of the radio jet (of order 10 7 yrs, ( 4.2 1, or 
even the duty cycle of jet activity of order 10 8 yrs estimated from 
observations of 'rejuvenated' radio sources (e.g., Schoenmakers 
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|et al.|[2000| . If the dissipation time is similar to the lifetime of 
the radio source, then the energy content of the gas will reflect 
the total energy injected by the radio source, including the early 
stages of radio activity when the jet was confined within the in- 
ner regions of the galaxy, and when the interaction was likely 
to be particularly efficient. If it is even longer, then it may in- 
fluence the interstellar medium and star-formation properties of 
the massive host galaxies of radio-loud AGN even on timescales 
beyond their 'active' phase. This is possible if the gas luminosity 
drops at the end of the jet activity period. In this case the dissi- 
pation timescales will be longer than our estimates. Addressing 
this question in detail would require observational constraints of 
the line emission in galaxies after their radio-loud phase (radio 
relics) and is beyond the scope of this paper. 

7. Molecular gas and star formation 

(Cold) molecular gas and star formation are closely related: gas 
surface densities and star-formation intensities (star-formation 
rate surface densities) show a close relationship over several or- 
ders of magnitude ('Schmidt-Kennicutt relation' Sch midt|1959| 
|Kennicutt|1998[ l. As we discussed previously, in H2 luminous ra- 
dio galaxies, however, star-formation rates appear low compared 
to the amount of available molecular gas. It is therefore inter- 
esting to investigate whether these galaxies fall onto the same 
star-formation law as 'ordinary' star-forming galaxies. 

7. 7. Constructing a Schmidt-Kennicutt-like diagram from 
PAH and CO line emission 

The best way of estimating star-formation rates is, generally 
speaking, the infrared continuum. However, in powerful radio 
galaxies the 24/im and 70/im flux scales with AGN power 
(Tadhunter et al. 2007), likely indicating that much of the dust is 
heated by the AGN rather than star formation. We will therefore 
rather use the PAH emission as an approximate tracer of star 
formation in our galaxies. For the radio galaxies of Ogl e et aL] 
(2010) we have direct measurements of the PAH fluxes from the 
IRS spectra, for star-forming galaxies we will use the 8/im lu- 
minosity instead, which in these galaxies is dominated by the 
7.7/zm PAH feature We use values corrected for the stellar con- 
tibution by Calzetti et al. ( 2007[ l; Roussel et al. ( 2007} and have 
already outlined in {3.1 how we translate 8jJ,m flux densities in 
PAH fluxes. 

Although PAHs are also not very robust tracers of the star- 
formation intensity generally, empirically they are found to be 
equally robust as the 24/zm flux when studying the integrated 
star formation rates in galaxies with high metallicities. (Calzetti 
|et al.|2007] ). The quality of our optical spectra of 3C 326 N is not 
sufficient to safely determine a metallicity from the absorption 
lines, but with its large stellar mass of order 10 M Q , 3C 326 N 
can safely be expected to have a high metallicity. The same ap- 
plies to H2-luminous radio galaxies generally. 

In Fig. [9] we show a Schmidt-Kennicutt-like diagram based 
on the PAH luminosity and CO surface brightness for 3C 326 N 
and other H2 luminous radio galaxies with published CO obser- 
vations ( |Ogle et al.|2010| lists H2 luminous radio galaxies with 
CO measurements in the literature), as well as for star forming 
and AGN host galaxies taken from SINGS ( |Roussel et al.|2007) . 
For galaxies where CO (and/or PAH) measurements are not spa- 
tially resolved, we assume a radius of 2.5 kpc each. We are using 
the same radius for PAH and CO emission, implying that molec- 
ular gas and PAH emission originate from the same region of the 
galaxy, which is astrophysically plausible. 



Using equations (2) and (6) o f|Calzetti et al.| ( [2007| ) we derive 
a calibration of star-formation intensity as a function of 8/im 
luminosity (vL v ) per kpc 2 or 



logT. SFR = -45.32 + 1.06 logSg 



>\lm 



(6) 



where T,sfr is the star-formation intensity in M yr~' kpc~ 2 
and S$n m the luminosity at 8/im in erg s _1 kpc -2 corrected for 
the stellar continuum (see Calzetti et al. 2007 , for details). We 
also use this relationship to re-calibrate the star-formation law 
of Kennicutt ( 1998) in terms of PAH luminosity, finding 



log%„, = 31.44+ 1.32 logE 



gas 



(7) 



where ~L gas is the molecular gas surface density in M 
kpc -2 . This relationship is somewhat less steep than the orig- 
inal |Kennicutt| law (which has a power-law index 1 .4), simply 
because PAH surface brightness and star-formation intensity do 
not scale linearly. 



7.2. Is the star-formation efficiency low in H2 luminous 
radio galaxies? 

The Schmidt-Kennicutt law is shown as the black solid line in 
Fig. [9] and is a good representation of the SINGS galaxies with 
CO emission line observations (Roussel et al. 2007). In strong 
contrast, 3C 326 N as well as the other H2 luminous radio galax- 
ies shows a pronounced offset by roughly 1-1.7 dex (a factor 10- 
50) towards lower PAH fluxes, but with a very similar slope. If 
the 'standard' CO conversion factor roughly applies, and if PAH 
emission traces star formation in a similar way as in 'ordinary' 
star-forming galaxies, then our result implies that the star for- 
mation in H2 luminous radio galaxies is about a factor 10-50x 
less efficient than in 'ordinary' star-forming galaxies. But even if 
these assumptions were not strictly appropriate, the basic result 
would persist, that the conditions of the interstellar medium in 
3C 326 N and other H2 luminous radio galaxies are significantly 
different from those in galaxies with 'ordinary' star formation 
properties. 

7.2.1. Caveats 

Fig.|9]is based on several assumptions, each of which introduces 
systematic uncertainties. This is why we prefer not to state a 
specific offset but rather give a range of possible offsets from the 
Schmidt-Kennicutt law. We will now discuss these uncertainties. 

One obvious caveat in our conclusion is that we only used 
the CO measurements to estimate a gas mass, neglecting the 
warm H2. Including the warm molecular gas in our mass esti- 
mate would only increase the offset of the H2 luminous galax- 
ies relative to the SINGS sample, for example by 0.3 dex for 
3C 326 N. Changes in the CO40-H2 conversion factor would 
have to be as large as the offset to move the H2 luminous galaxies 
onto the nearby relationship, and would again not change the ba- 
sic conclusion, namely that the physical gas conditions in these 
galaxies are markedly different from those in 'ordinary' star- 
forming galaxies. In this case, the CO-to-H2 conversion factor 
would have to be lower by a factor 30 than in the SINGS galax- 
ies, and a factor 6 lower than in ULIRGs ( |Downes & Solo mon 
[T998T >. 

Another possible caveat is that we adopted a fiducial radius 
of 2.5 kpc for the galaxies where we only have integrated mea- 
surements (empty red triangles in Figure |5J. Varying this radius 
by large amounts (greater than factors 2-3) will mildly affect 
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the position of individual galaxies, but is not sufficient to put 
them onto the 'ordinary' Schmidt-Kennicutt ridge. This limita- 
tion can easily be overcome by collecting larger samples with 
high-resolution, high-quality CO observations. Galaxies shown 
as red filled triangles have spatially-resolved CO measurements 
and are not affected by this effect. We adopted the same radii 
for PAH and CO emission, which corresponds to the astrophys- 
ical assumption that both are related to the same star-forming 
regions. 

We may also wonder whether this offset could reflect a net 
deficit of PAH emission? We will rely on the close analogy with 
the H2-luminous shock in Stephan's Quintet to address this con- 
cern. For Stephan's Quintet, [Guillard et al. ( 2010) > show that the 
PAH and mid-IR continuum associated with the warm H2 emis- 
sion is consistent with the expected emission from the warm 
molecular gas, assuming that the dust-to-gas mass ratio, and the 
fraction in PAHs and very small grains have roughly the same 
values as in the Milky Way. 

The ratios between H2, CO and PAH luminosities in 3C 326 
are very similar to those measured in Stefan's Quintet (Figure[6]i, 
so that similar arguments may apply to 3C 326 N. For a stan- 
dard (Galactic) dust mass fraction in PAHs, the weakness of the 
PAH emission implies a UV field comparable to that observed in 
the intergalactic shock in Stefan's Quintet. The UV intensity in 
Stephan's Quintet, measured with GALEX, is not very different 
from that in the Solar Neighborhood (1.4 in Habing units^ji. We 
used Starburst99 (Leithe rer et al.|19 99 ) to deduce an upper limit 
on the UV radiation field from the measured up per limit on the 
star-formation rate (0.07 Mq yr _1 , inferred by |Ogle et al.|2 007 
from the PAH and mid-infrared dust continuum of 3C 326 N). 
Assuming that the molecular gas is within a sphere of 2.5 kpc 
radius (the approximate extent of the CO emission), we find an 
upper limit on the UV radiation field of about 1 Habing unit. 
Thus, the weakness of the PAH and dust mid-IR emission is con- 
sistent with the low star-formation rate, and does not require the 
molecular gas to be PAH poor. Measuring the dust mass and tem- 
perature directly from the far-infrared luminosity is now possible 
with Herschel. This will allow for more robust estimates of the 
star-formation rates, and give direct constraints on the CO-to-H2 
conversion factor. These measurements will be critical to better 
estimate the star-formation efficiencies in H2 luminous galaxies. 



fragmentation and gas collapse and ultimately star formation. 
In Cen A (NGC 5128) |Neumayer et al.| ( |2007] > find that the ro- 
vibrational H2 emission lines are broader than expected from the 
stellar velocity dispersion, although they have an overall sim- 
ilar velocity field. This would certainly agree with a scenario 
where a small fraction of the energy injected by the AGN into 
the multiphase interstellar medium of the host galaxy will finally 
contribute to increasing the turbulence in the molecular disk. 
Detailed follow-up observations of a sufficiently large sample 
of H2-luminous galaxies are certainly necessary to substantiate 
this speculation. 

We did not find evidence for positive AGN feedback en- 
hancing the star formation efficiency relative to 'ordinary' star 
forming galaxies, as propos ed by, e.g.,|van Breugel et al. ( 1985|l; 
Begelman ~& Cioffi|p9g9| l; |Mellema et al.|(|2002);|Fragile et al. 
P004);|Croft et al.|(|2006|;|Silk & Norman[(|2009|l.|van Breugel 
l et al.| (|1985|);|Schirninovicn et al.| ( |1994) ; |Croft et al.| (|2006 ); 
Elbaz et al. (2009) claim evidence for star formation associated 
with radio jets outside the host galaxy in a few individual nearby 
systems, which may be triggered by the jet. Indeed, turbulence 
can trigger star formation locally where gravitationally-bound 



giant molecular clouds are otherwise not able to form ( Klessen 
et al. |20()0 1. If positive feedback dominated in our galaxies, we 



would expect an offset to higher star-formation intensities for a 
given gas surface density in Fig. [9| which is not the case, as all 
radio galaxies are shifted towards lower star-formation efficien- 
cies. 

Our overall results suggest that the suppression or enhance- 
ment of star formation depends critically on the details of the 
'micro-'physics of the gas on small scales, which is beyond the 
capabilities of current numerical models and requires detailed 
observations of the multiphase gas in AGN host galaxies, which 
are becoming possible only now with Spitzer, Herschel, ALMA, 
and JWST In particular, it is unclear whether our results can 
be easily extrapolated to (radio) galaxies at high redshift, which 
have copious amounts of molecular gas traced through CO fine 
emission (e.g., Papadop oulos et al.|2000||De Breuck et aL|2 003 



2005 [Nesvadba et al. 2009), but whose gas conditions are likely 



very different from those in 3C 326 N. If the kinetic energy in 
these galaxies can be dissipated rapidly enough to allow for gas 
collapse and star formation, then feedback may be positive. 



7.2.2. Implications 

We are not the first to note that star formation in radio galaxies 
with CO detections appears suppressed relative to radio-quiet 
galaxies. Oku da et al.| ( |2005) and Papadopoulos et al. ( 2007) 
found the same effect for 3C31 and 3C293, respectively. Koda 
|et al.| ( |2005] l identify star-forming knots embedded in the molec- 
ular disk of 3C31, pointing out that, although the disk appears 
globally stable against gas collapse and star formation, this is 
not sufficient to quench all star formation. The position of 3C3 1 
in Fig. [9] suggests that, nonetheless, the overall star formation 
efficiency in this galaxy is considerably lower than in 'ordinary' 
star-forming galaxies. 

Our previous discussion of how the radio source influences 
the energy budget of the interstellar medium of 3C 326 N (Sj6]l 
suggests that heating through the radio source may be an attrac- 
tive non-gravitational mechanism, at least in galaxies that host 
radio-loud AGN, which may enhance the turbulence and heat- 
ing in a molecular disk, preventing or suppressing gravitational 

2 The Habing field is equal to 2.3 ergs~'cm~ 2 at A = 1530 A 
(Habingl 1968b 



8. Implications of this scenario for galaxy 
evolution 

The energy injected by powerful AGN into the interstellar 
medium and halo of galaxies may play an important role in de- 
termining the characteristics of galaxies in the local Universe. 
For example, the high metallicities, luminosity-weighted stellar 
ages, and relative abundances of a elements relative to iron sug- 
gest rapid, and truncated star formation in the early Universe 
followed by a long phase of (mostly) passive evolution (e.g., 
Pipino & Matteucci 2004) for massive galaxies. These obser- 
vations are consistent with a phase of powerful AGN driven 
outflows in the early Universe. In fact, studies of powerful ra- 
dio galaxies at 2 with rest-frame optical integral-field spec- 
troscopy reveal energetic outflows which have the potential of 
removing a significant fraction of the interstellar medium of 
a massive, gas-rich galaxy within the short timescales neces- 
sary to explain the super-solar [a/Fe] ratios and old luminosity 
weighted ages observed in massive ellipticals ( |Nesva dba et al. 
2006 , 2007] |2008[ ). This "quenching" of star formation by pow- 
erful radio jets is in broad agreement with theoretical models 
explaining the characteristics of massive galaxies. However, in 
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addition to this quenching phase, the radio-loud AGN may also 
assist in maintaining low star-formation rates over cosmological 
timescales by inhibiting subsequent gas cooling - "maintenance 
phase". This last phase is important since, even in the absence 
of mergers, subsequent gas infall and return from the evolving 
stellar population will replenish the reservoir of cold gas which 
must be prevented (at least partially) from forming stars. 

In the scenario we presented above, roughly equal amounts 
of energy are being dissipated through turbulence, and are pow- 
ering an outflow of warm neutral and ionized gas. Gas that does 
not escape from the halo of the host galaxy, will likely cool 
and rain back onto the galaxy. For cool-core clusters, Revaz 



formation rate of SFR^jL syn 
with the current SFR^ 3 



^4 M yr 1 is broadly consistent 

--1 



et al. (2008); Salome et al.|(2008) propose that molecular fila- 



ments may represent gas that has been lifted from the galaxy, has 
cooled, and may be "raining back down" into the deepest part of 
the gravitational potential. A similar mechanism may apply to 
individual galaxies like 3C 326 N. In addition, the stellar popu- 
lation will continue to feed the interstellar medium with material 
which will also cool and descend into the potential well, if the 
material will not have been ablated previously. This may once 
again fuel an AGN, leading to a self-regulating AGN fueling cy- 
cle, even in the absence of galaxy merging. 

We have so far focused our discussion on 3C 326 N which is 
the most extreme H2 luminous radio galaxy known and there- 
fore well suited to study the physical mechanism responsible 
for the H2 line emission. However, galaxies are complex astro- 
physical objects, where a wide spectrum of physical and astro- 
physical phenomena may act in parallel. In order to illustrate 
how different formation histories may lead to different obser- 
vational signatures within a common physical framework, we 
will in the following contrast 3C 326 N with 3C293 at z=0.045. 
We have chosen 3C293 because it has particularly luminous CO 
line emission (Evans et al.|1999| and unlike 3C 326) and with 
CO line ratios suggesting the molecular gas is highly excited 
( Papadopo ulos et aL|2008). It als o has a significant outflow of 
neutral gas ( [Morganti et ai|2003H2005||Emonts et al.|2005) . Its 
H2 emission-line luminosity, warm molecular gas mass and op- 



tical line ratios are very similar to those of 3C 326 N (Ogle et al. 
2010} , which suggests that 3C293 may be a 'gas-rich analog' 
and thus a foil of 3C 326 N. 



8. 1. Differing Star-formation Histories: Clues to the 
importance of molecular gas content? 

3C 326 N and 3C293 have very different star-formation histo- 
ries. The spectrum of 3C 326 N is consistent with an old stellar 
population with an age of > 10 10 yrs. At most a small fraction of 
the total stellar mass, < 5%, formed recently (~ l-2x 10 9 yrs; 
Fig.|7j. Within the 3" aperture of the SDSS fiber, this would cor- 
respond to 3 x 1O 9 M of stars formed during the last 1-2 Gyrs, 
implying a star-formation rate of about 1M yr~ 1 for a constant 
star-formation history. We note that this is more than a factor 10 
greater than the upper limit on the current star-formation rate of 
<0.07 M0 yr -1 estimated from the infrared emission. 

For 3C293 we find a significantly more complex star- 
formation history (see also Tadhunter et al. 2005 1. Namely ~ 
80% of the stellar mass of 3C293 was formed at high redshift, 
20% were formed in a more recent star-forma tion episode 

found sim- 



about 1 - 2 x 10 9 yrs ago (jTadhunter et al. 



2005 



ilar results). In addition, the equivalent widths and ratio of the 
Ca H+K absorption lines require a very young population with 
an age of ~ 10 7 yrs, consistent with the starburst implied by 
far-infrared luminosity (Papadopoulos et al. 2008). The current 
burst contributes <1% to the total mass of 3C293, and the star- 



7 M Q yr estimated from the infrared 
luminosity of 3C293 (Papadopou los et al.||2008| l. To evaluate 
whether this fit is unique, we tested alternative star-formation 
histories with only one or two of the three populations. Fig. [7] 
illustrates that the three-component fit indeed is a significantly 
better representation of the data. The population synthesis fits 
also yield stellar mass estimates, but we need to correct for the 3" 
fibre size of the SDSS, which is much smaller than the size of the 
galaxies. Doing this, we find stellar masses of M^^L = 3 x 10 11 

M and M®% lar = 7 x 10 1 'M for 3C 326 N and 3C293, respec- 
tively. 



8.2. Origin of the gas 

3C293 has a ~ 10 x larger reservoir of gas traced through CO 
line emission than 3C 326 N. What is the likely origin of the 
gas in these galaxies? Both galaxies have nearby galaxies with 
which they are perhaps interacting. For the pair 3C 326 N/S, we 
find a very low gas mass and red optical colors for both galaxies, 
suggesting this is a gas-poor system. Indeed our estimate of the 
stellar and gas mass suggests a gas fraction of only ~ 1%. The 
star-formation history of 3C 326 N appears quiescent and does 
not indicate that the galaxy was undergoing a starburst due to the 
interaction. Nonetheless, we may expect the accumulation of a 
significant amount of cold gas due to mass loss from its stellar 
population (~1 Mq yr -1 ) as well as perhaps accretion from the 
surrounding halo. How much gas will this likely be? 

1994) ) to estimate 



We use the analysis of |Kennicutt et al. 



the amount of gas returned by a single-age stellar population 
into the interstellar medium over 10 Gyr of passive evolution. 
Kennicutt et al. ( 1994 1 find that depending on the details of the 
star-formation history and initial mass function, as much as 30- 
50% of the stellar mass formed in the burst may be returned into 
the ISM during 10 Gyr. From their Fig. 8 we estimate that the 
gas return after the first 0.5-1 Gyrs will be of order a few per- 
cent of the stellar mass. For the stellar mass of 3C 326 N, and 
ignoring any other mechanism like subsequent star formation or 
feedback, this would correspond to a total gas mass of several 
10 9 M up to~ 1 x 10 10 M , factors of a few more than the 
amount of cold and warm gas we observe in 3C 326 N. In addi- 
tion, with an infall rate of ~ 0.1-1 M yr -1 of gas accreted from 
the halo (which appear s plausible for early and late-type galax- 
ies; Sancisi et al.|2008 1 over 10 10 yrs, this gas mass may roughly 
double. An infall rate of ~ O.1M yr -1 of neutral gas has also 



been estimated by Morganti et al. ( 2009 1 from the observation of 
a narrow, redshifted HI cloud seen in absorption and emission in 
the nearby, radio-loud early-type galaxy NGC315. 

Hence, each of these processes alone could easily explain the 
observed gas mass of 3C 326 N. For 3C293 however, similar es- 
timates of stellar mass loss and gas accretion would be inefficient 
to explain the amount of gas necessary for forming the ~ 10 11 
M of stellar mass during the last Gyr, which is certainly consis- 
tent with the assumption of a merger-triggered starburst for this 
galaxy. So the difference could be in the "mass accumulation 
histories" of these two AGN. 

8.3. Maintaining the low fraction of recent star-formation 

Despite these plausibly different accretion histories, both of the 
galaxies apparently have low star-formation efficiencies. We es- 
timated that a total of 1 — 2 x 10 10 M of warm/cold gas would 
accumulate in a massive early-type galaxy such as 3C 326 N 
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over a Hubble time while 3C293 acquired its gas during an in- 
teraction/merger. Is it plausible that the gas in both sources will 
be heated and substantial fractions removed by the mechanical 
ene rgy of the jets? 

|Best et aL] ( |2005| ) find that the number of radio sources is a 
strong function of radio power and stellar mass of the host galaxy 
for a sample of 2000 early-type galaxies with FIRST and NVSS 
observations taken from the SDSS. For example, for galaxies 
with stellar mass of ~ 3 x 10 11 M Q they estimate that ~ 10% of 
all galaxies have radio powers > 10 24 W Hz ', and propose that 
this may represent a duty cycle, where any given galaxy of this 
stellar mass will host a similarly powerful radio source for 10% 
of the time. 

Statistically speaking, over 10 Gyrs, a given AGN in a few 
x 10 11 M Q galaxy would be active for a total of ~ 10 9 yrs (this 
total 'activity period' may consist of many radio-loud phases 
with i nter mediate, quiescent phases). With the same reasoning 



as in 



4.2 



we estimate that a power > 10 24 W Hz -1 measured at 
1 .4 GHz will correspond to a kinetic power of order 10 43 erg s~ 1 , 
equivalent to a total energy of a few x 10 59 erg for a total activity 
time of 10 9 yrs (about 10% of a Hubble time). A few x 10 58 erg 
of energy are necessary to unbind 10 10 M in gas fro m the po- 
tential of a galaxy with few x 10 M in stellar mass ( jMorganti 
et al.|2005 1 [Nesvadba et al.|2006||20"08] i. Based on these simple 
energy considerations, it appears overall plausible that AGN in 
massive galaxies may regulate gas cooling through the energy 
injected by their radio sources over a Hubble time as expected 
for the "maintenance" mode. 

However, in addition to the energy requirement, we also have 
a timescale requirement, since gas cooling must be inhibited in- 
between phases of radio-loud AGN activity, if AGN are to shape 
the properties of the ensemble of (massive) galaxies. Several ar- 
guments may suggest that this could in deed be possible. First, 



the dissipation timescales we found in S6.3 of order 10 yrs 



imply that significant fractions of the molecular gas will remain 
warm during all of the lifetime of the radio jet, and possibly also 
for a significant time afterwards (if the upper estimate, albeit 
not very precise, is more appropriate). Interestingly, studies of 
'rejuvenated' jet activity suggest timescales of about 10 8 years 
between activity episodes ( Schoenmakers et al. 2000), similar to 
the upper range of the H2 dissipation time. Second, the mechani- 
cal energy injected by radio-loud AGN is sufficient to accelerate 
and perhaps remove significant amounts of gas. Even if parts of 
the entrained material will not reach escape velocity, its density 
will decrease due to the larger volume. The mean free path of 
particles between collisions will be greatly increased, making 
the dissipation timescales accordingly longer and delaying the 
time until this gas will cool and rain back onto the galaxy. Third, 
the expelled gas may be reheated by mixing with the hot halo 
gas. Fourth, the molecular gas we observe may only be tran- 
sitory, where the molecules are formed by the ram pressure of 
diffuse gas (Glover & Mac Low 2007) as a result of the dissi- 
pation of the turbulent energy. This process does not necessarily 
lead to the formation of gravitationally bound clouds, in which 
case the H2 destruction through photodissotiation may dominate 
over the H2 formation, when the turbulent and thermal pressure 
in the cocoon drop. 

9. Summary 

We presented a detailed analysis of the physical gas conditions 
in the nearby powerful, H2-luminous radio galaxy 3C 326 N, 
which does not show the signatures of active star formation in 
spite of few x 10 9 M Q in molecular gas. Our main goals were to 



investigate the gas conditions in this galaxy, to unravel the astro- 
physical mechanism which powers the remarkably bright mid- 
infrared H2 emission, and to elucidate why this galaxy is not 
forming stars in spite of a considerable reservoir of molecular 
gas. To this end, we combined newly obtained CO(1-0) millime- 
ter spectroscopy of 3C 326 N obtained at the Plateau de Bure 
Interferometer with the Spitzer-IRS spectroscopy by Ogle et al. 
(2007] |20T0| and optical spectra from the SDSS. Our main con- 
clusions are as follows: 

(1) We compare gas masses for warm and cold gas, estimated 
from Spitzer mid-infrared spectroscopy and IRAM millimeter 
imaging spectroscopy. We find that most of the molecular gas in 
3C 326 N is warm (~2x 10 9 M , compared to 1.5 x 10 9 M 
estimated from our CO(1-0) observations for a "standard" CO- 
to-H2 conversion factor). This ratio of warm to cold gas mass is 
about 1 to 2 orders of magnitude larger than that found in star- 
forming galaxies, indicating that the gas is in a distinct physical 
state, which may account for the low star-formation efficiency. 

(2) We introduce a new 'molecular' diagnostic diagram based on 
the pure-rotational H2, CO and PAH emission to show that most 
of the line emission from molecular gas in 3C 326 N is not pro- 
duced by UV heating. We argue that the gas is likely to be pow- 
ered by the dissipation of mechanical energy through shocks. 
The same is found for the ionized gas. 

(3) Interstellar Na D absorption marks an outflow of neutral gas 
at velocities of up ~ —1800 km s _1 , and most likely mass and 
energy loss rates of 30—40 M yr _1 and ~ 10 43 erg s _1 , respec- 
tively. These values are similar to those found in starburst-driven 
winds, but this is the first time that such an outflow is detected 
in the Na D line of a galaxy which does not have strong star 
formation. The star-formation intensity in 3C 326 N is orders of 
magnitudes below what is necessary to drive a wind. Similarly, 
the line profile cannot be explained through the interaction be- 
tween 3C 326 N and 3C 326 S. 

(4) Based on these observations we propose a scenario where 
the outflow and H2 line emission are intricately related. It repre- 
sents an extension of the well-explored 'cocoon' -model of inter- 
actions between radio jets and the ambient gas, and includes the 
physics of the molecular gas. In this scenario,the H2 line emis- 
sion is powered by turbulence induced within dense clouds that 
are embedded in the expanding cocoon, with a dissipation time 
of ~ 10 7 ~ 8 yrs. Dissipation times of 10 8 yrs are in rough agree- 
ment with the duty cycle of jet activity estimated from rejuve- 
nated radio sources. 

(5) Comparing PAH and CO surface brightness (in analogy to 
the "Schmidt-Kennicutt" diagram), we find a significant off- 
set between 3C 326 N and other H2-luminous galaxies towards 
lower PAH surface brightnesses. This may suggest that star- 
formation efficiencies in these galaxies are lower by roughly a 
factor ~ 10 — 50 than those in 'ordinary' star-forming galaxies. 

(6) Generalizing our results for 3C 326 N we find that outflows 
during similar radio-loud episodes in massive galaxies may bal- 
ance the secular supply in cold gas through accretion and the 
mass return from evolved stars over a Hubble time. If radio- 
activity is a common, but episodic property of most massive 
early-type galaxies as suggested by |Best et al. ( 2006[ ), then the 
long dissipation times of up to ~ 10 8 yrs may have consequences 
for the population of massive galaxies as a whole. 

We emphasize that this analysis can only be the first step 
towards an understanding of how radio-loud AGN regulate the 
physical conditions of the molecular gas and ultimately star for- 
mation in massive galaxies. It illustrates the importance of fol- 
lowing a multi-wavelength approach if we seek to relate the 
physics and astrophysics of AGN feedback with the physical 
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conditions of the interstellar medium including molecular gas 
in the host galaxy and the regulation of star formation. 
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Fig. 1. Continuum-free CO(l-O) emission-line morphology of 3C 326 N (white contours) superimposed on the SDSS R-band image 
of 3C 326 N and 3C 326 S. Contours are given for 3, 4, and 5(7. The inset shows the same SDSS R-band with the line-free 3 mm 
continuum morphology superimposed as white contours. We used the continuum morphologies to align the millimeter and optical 
data, assuming that the AGN reside in the nuclei of the galaxies. With this alignment, we find an offset of 0.7" between the CO line 
emission and continuum (millimeter and optical) which, given the faintness of the line emission and much greater beam size, is not 
significant. The ellipse in the lower left corner shows the beam size and position angle in both images. 
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Fig. 2. CO(1-0) millimeter spectrum of 3C 326N (red solid line), and H/3 line profile (blue solid line). The red dashed line shows 
our best fit to the CO(1-0) line profile. Both spectra are continuum-subtracted. 
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3C326 H 2 Excitation Diagram (n H =10 3 cm -3 ) 
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Fig. 3. Observed H2 excitation diagrams for 3C326 (1<7 error bars) and model results of magnetic shocks (see Guillard et al.|2009 



for details), for pre-shock densities of n# = 10 3 (upper panel) and 10 4 cm 3 (lower panel). Dashed lines show the contribution of 
each C-shock velocity to the best-fit model. The red line shows the combination of these 3 shocks. 
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Fig. 4. Multiple-component fits to the [OIII]AA4959,5007 and [OI]AA6300,6363 line emission in 3C 326 N. The broad component 
is best seen in the [OIII]A5007 emission line. We assumed the same redshifts and line widths for all lines. 
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Velocity [km 3''] rest wavelength [A] 

Fig. 5. left: NaD absorption line in 3C 326 N. In the upper panel, the spectrum is shown in black, and the red line shows the stellar 
continuum for our best stellar population fit. The two NaD components are not resolved at R=1800. In the mid panel, we show the 
spectrum with the best-fit stellar continuum subtracted. The green line in the mid panel marks our absorption line fit with HelA5876, 
and two Gaussian NaD components, one at the systemic velocity, another one with an offset of -350 km s representing the blue 
wing. The blue line shows the Hel line and the stellar NaD component. In the lower panel, we show the spectrum with the stellar 
continuum and the Hel emission-line subtracted (upper spectrum) and the residual after subtraction of Hel, the stellar component, 
and the NaD absorption line corresponding to our best-fit wind model. The residual spectrum is shifted by an arbitrary amount along 
the ordinate. The abscissa gives the velocity offsets relative to a systemic redshift of z=0.090. The Na D profile shows evidence for 
absorption with blueshifts of up to ~ —1800 km s _1 , and potentially, at lower signal-to-noise levels, of up to —6000 km s . right 
The same for Mgb, a line which does not have a strong interstellar component. The residual (lower panel) is consistent with noise 
after subtraction of the stellar component and the NIA5199 line. 
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Fig. 6. Molecular diagnostic diagram to distinguish between AGN and star formation based on PAH bands, CO(1-0) and mid- 
infrared H2 line emission. Lines mark PDR models with different assumptions for the UV radiation fields. Star-forming galaxies 
from the SINGS survey fall within the portion of the diagram spanned by PDR models (black lines), unlike Stephan's Quintet, 
3C 326 N and a number of other H-2 luminous radio galaxies with published CO(1-0) observations. Consequently, the H2 line 
emission in these galaxies cannot be powered by energetic photons produced in star-forming regions. Values for the SINGS galaxies 
are taken from Roussel et al. ( 2007| l. H2 and PAH fluxes for the H2 luminous radio galaxies are taken from Ogle et al. (2010), who 
show the relationship between H2-to-PAH 7.7/im versus 24/im continuum luminosity. 
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Fig. 7. Wavelengths near the 4000 A break are most sensitive to the star-formation history of 3C 326 N (upper panel) and 3C293 
(lower panel). The black line marks the SPS S spectrum in each panel. Red, blue, and green lines mark the Starlight spectral fits 



for different star-formation histories, see £8.1 for details. Upper Panel: Red - Stellar population older than 10 10 yr. Blue - Old and 
intermediate-age (few 10 9 yrs) stellar population. Both scenarios fit the data equally well. Lower Panel: Red - Stellar population 
older than 10 yr. Green - Old and intermediate-age (few 10 8 ~ 9 yr) population. Blue - Old , intermediate-age , and young (10 7 
yr) stellar population. A single star-formation episode 10 10 yrs ago does not appear as a good fit to the data; adding a young stellar 
population (consistent with the observed starburst) improves the fit to the Ca H+K lines, although significant differences remain, 
especially for the blue component. 
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Fig. 8. This diagram outlines the energy flow within the hot "cocoon" that we propose to account for the HII and H2 line emission 
and the wind. The flow is powered by the interaction of the jet with the multiphase interstellar medium of the host galaxy (arrow to 
the left). It feeds three energy reservoirs: (boxes from left to right:) the thermal energy of the hot plasma surrounding the relativistic 
radio source, the turbulent kinetic energy of the clouds embedded in the plasma and the thermal energy of the warm HII and 
H2 gas. The arrows at the bottom represent the "cocoon" energy loss by radiative cooling, that to the top the transformation of 
thermal energy into bulk kinetic energy of the outflowing gas. The dissipation of the turbulent kinetic energy powers the H2 and 
HII line emission. The high pressure plasma expands on a time scale shorter than its cooling time through X-ray emission driving 
a multiphase outflow. Our observational results suggest that the kinetic power of the outflow and the power radiated by the HII and 
H2 gas are both ~lx 10 43 erg s . The power input represents ~ 10% of the mechanical power of the jet. 
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Fig. 9. PAH intensity as a function of molecular gas mass surface density measured from the CO(l-O) emission-line intensity. The 
black solid line shows the Schmidt-Kennicutt law (pCennicutt 1998 ), where the star-formation intensity is translated into a PAH 
intensity using the results of |Calzetti et al.| p007 ). SINGS star-forming and AGN host galaxies from Roussel et aL]( |2007[ ) are shown 
as small black triangles and small blue circles, respectively, and fall onto the relationship. In stark contrast, 3C 326 N and other H2 
luminous radio galaxies (big red triangles; PAH fluxes are taken from |Ogle et al.pOlO l have a marked offset by roughly factors 
10-50 towards lower PAH intensities for a given molecular gas mass (black hatched lines), but are consistent with having a similar 
slope as the 'standard' Schmidt-Kennicutt law. This offset would become larger if we had included the warm molecular gas (e.g., 
by 0.3 dex for 3C 326 N). (Filled red triangles mark galaxies with spatially resolved CO detections, empty triangles mark galaxies 
with integrated measurements. Arrows mark galaxies with sensitive upper limits.) 
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Aperture 


v obs 


Peak 


Width 




[GHz] 


[mJy] 


[km/s] 


2.5" x 2.1" 


105.744 


1.1 ±0.2 


351±98 


5"x5" 


105.744 


1.5±0.3 


351±98 



Table 1. Observational results: CO(l-O) line-fit parameters of (1) a spectrum integrated over a region of 2.5" by 2.1" and (2), a 
spectrum integrated over a 5"x5" aperture corresponding to the slit width of the IRS spectrum (dashed line in Fig. [2} 



Aperture 


redshift 




U 


Fcont 






[Jy kms -1 ] 


[Kkms- 1 pc 2 ] 


[mJy] 


2.5" x 2.1" 


0.0901 


0.5±0.1 


2±lxl0 s 


0.99 


5"x5" 


0.0901 


1.0±0.2 


3.8±lxl0 8 


0.99 



Table 2. Observational results: Intensity of the CO(1-0) emission line after subtraction of a point-like continumm source with a 
fixed value of 0.99 mJy at the position of the cm radio-source (Rawlings et al. 1990). L' co was computed with the formula of 
Solomon et al. (1997). Results are given for (1) a spectrum integrated over a region of 2.5" by 2.1" and (2) a spectrum integrated 
over a 5"x5" aperture corresponding to the slit width of the IRS spectrum (dashed line in Fig. El. 



line 


^rest 


redshift 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


[Oil] 


3727.00 


0.09019 


623 


9.62 


Hj8 


4861.30 


0.09017 


610 


2.44 


[om] 


4959.90 


0.08991 


609 


0.58 


[om] 


5006.90 


0.09011 


609 


1.77 


[OI] 


6300.30 


0.09013 


612 


2.22 


[OI] 


6363.80 


0.09013 


612 


0.53 


[Nil] 


6548.10 


0.09014 


601 


2.38 


Ha 


6563.30 


0.09007 


596 


9.30 


[Nil] 


6584.10 


0.09003 


601 


7.17 


[SII] 


6716.40 


0.09019 


605 


4.86 


[sn] 


6730.80 


0.09019 


606 


3.68 



Table 3. Emission-line fluxes in 3C 326 N. Column (1) - Line ID. Column (2) - Rest-frame wavelength in A. Column (3) - Redshift. 
Column (4) - Full width at half maximum in km s , corrected for an instrumental resolution of 167 km s , corresponding to a 
spectral resolving power of R=1800. Uncertainties are a few 10s of km s . Column (5) - Integrated line flux in units of 10~ 15 erg 
s _1 cm~ 2 . All values are derived after subtracting the best-fit stellar continuum (see text for details). Uncertainties are a few % for 
all lines. For [OII]AA3727, additional uncertainties are due to the unknown line ratio between the two components of the doublets. 
These are insignificant for the FWHM and flux estimate, and amount to ~ 70 km s for the redshift, which however is not one of 
the important quantities for our analysis. 



Line ratio 


observed 


shock 


shock & precursor 


(1) 


(2) 


(3) 


(4) 


HjS/Ha 


0.26 


0.33 


0.33 


[OIII]/H/3 


0.72 


0.82 


1.48 


[OI]/Hce 


0.24 


0.22 


0.13 


[NII]A6584/Ha 


0.77 


0.68 


0.47 


[Sn]/Ha 


0.61 


0.52 


0.32 


[SII]A16716,6731 


1.32 


1.23 


1.24 


[NeII](12.81,um)/Ha 


0.46 


0.45 


0.28 


[Nell ] ( 1 2 . 8 1 ^m)/[NeIII] ( 1 5 .6/im) 


>1.8 


2.2 


0.73 



Table 4. Measured emission -line ratios in 3C 326 N and expected ratios for a pure-shock model and a model assuming a shock and 
precursor ( Allen et al.|2008| . The line ratios are consistent with shock velocities of v=250 km s _1 , a pre-shock density of 1 cm" 1 , 
and a magnetic parameter of ~&lyj\n) < 1 /iG cm~ 3 / 2 . Column (1) - Line ID. Column (2) - Observed line r atio. Column (3) - Pure 
shock model. Column (4) - Shock and precursor. [Nell] and [Nelll] fluxes are taken from Ogle et al. (2010 1. 
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Molecular Gas 
( l ) 


Ionized Gas 
\ L ) 


Radio source 
\->) 


ACjJN A-ray 


obs. luminosity 


41.9 


41.3 


44.6 


40.6 


bol. luminosity 


42.0 


43 






Kinetic energy 


57.5 


55.7 


59.8 




Mass 


9.5 


7.3 






~diss 


8.2 


7.2 







Table 5. Overview of the luminosities of molecular line emission and Ha (in erg s ), kinetic energies (in erg s ), masses (in M ), 
and dissipation times (in yrs). We give the logarithms for all values. Column (1) - Molecular gas measured with IRAM and Spitzer. 
The dissipation time is given under the assumption that all of the kinetic energy is dissipation through molecular line emission. 
Column (2) - Warm ionized gas (the observed luminosity is that of Ha). The dissipation time is given under the assumption that all 
of the UV/optical line emission is due to the dissipation of kinetic energy. Column (3) - Radio luminosity measured at 327 MHz. 
The kinetic energy was derived by Willis & Strom (1978). Column (4) - X-ray luminosity of the nuclear point source given by Ogle 
|etaL| ( |20T0) . 



Table 6. MHD shock model parameters and predicted H2 line fluxes" for 3C 326. 



r, b 


Vs c 








H 2 Line Fluxes [10 _lii W mr 1 } 










[cm" 3 ] 


[kms- 1 ] 


S(0) 


S(D 


S(2) 


S(3) 


S(4) 


S(5) 


S(6) 


S(7) 




4 


0.74 


4.51 


0.18 


6.1 x 10~ 3 


3xl0~ 4 


8.6xl0~ 4 


2.9 xlO -4 


7.5 xlO -4 


5.43 


13.02 


10 4 


15 


0.09 


3.22 


3.19 


10.36 


1.35 


0.68 


0.02 


2.2xl0~ 3 


18.92 


20.55 


34 


4.2xl0" 3 


0.19 


0.29 


2.32 


1.38 


5.33 


1.58 


2.99 


14.09 


16.35 




Sum 


0.84 


7.92 


3.66 


12.69 


2.73 


6.01 


1.60 


2.99 


38.44 


49.92 




4 


1.61 


3.74 


0.018 


1.3 xlO" 3 


4.3 xlO" 4 


1.2 xlO" 3 


4xl0~ 4 


10- 3 


5.37 


10.6 


10 3 


24 


0.08 


3.01 


2.44 


5.09 


0.46 


0.17 


7x 10~ 3 


6.6xl0~ 3 


11.3 


11.6 


40 


0.02 


1.01 


1.32 


7.57 


2.73 


5.24 


0.68 


0.56 


19.13 


20.07 




Sum 


1.71 


7.75 


3.78 


12.66 


3.2 


5.42 


0.68 


0.56 


35.77 


42.27 


Spitzer Obs. 


2.0±0.3 


7.5±0.5 


4.1±0.4 


12.6±0.4 


3.1±0.6 


5.0±2.5 


2.5±1.0 


2.9±1.0 


39.7±2.3 





° This table lists the shock model velocities, the predicted H2 rotational line fluxes, and bolometric luminosities associated with each shock 
velocity components. Measured H2 fluxes are taken from Ogle et al. |2007 2010 1. 
Preshock hydrogen density. 
MHD Shock velocity. 

d Sum of the H 2 S(0) to S(7) rotational lines in 10~ 18 W irT 2 
e Sum over all the lines (bolometric luminosity of the shock) 



nH 


V s 


Mass flow 


tcool 


M eas 


(1) 


(2) 


(3) 


(4) 


(5) 


10 4 


4 


70100 


19050 


1.3 x 10 y 


10 4 


15 


5950 


3320 


2.0 x 10 7 


10 4 


34 


730 


1000 


7.3 x 10 5 


10 3 


4 


41050 


65000 


2.7 x 10 y 


10 3 


24 


1100 


11500 


1.3 x 10 7 


10 3 


40 


630 


5400 


3.4 x 10 6 



Table 7. MHD shock model parameters, mass flows and cooling times. Column (1) - Preshock density in cm 3 . Column (2) - 
Shock velocity in km s . Column (3) - Mass flow in M yr~'. Column (4) - Cooling time in yrs. Column (5) - Gas mass in M Q . 



